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A  differential  reflectometor,  which  is  capable  of  measuring 
small  differences  in  reflectivity  between  two  samples,  is  described. 
Among  the  possible  applications  of  this  device  are:   to  reveal  the 
structure  in  the  spectral  reflectivity  of  metals  and  alloys;  to  study 
changes  in  reflectivity  caused  by  surface  reactions  and  environmental 
changes,  and  to  investigate  metallurgical  reactions. 

The  instrument  has  been  used  in  the  study  of  compositional 
modulation  of  copper-based  alloys  at  room  temperature  in  the  spectral 
range  between  1.8  and  6.2  eV. 

For  copper-nlck«l  rilloys,  a  decrease  in  reflectivity  in  the 
low-energy  region  was  observed  with  increasing  concentrations  of  nickel. 
No  significant  shift  of  the  fundamental  edge  at  2.1  eV  in  pure  copper 
was  detected  when  nickel  was  added  to  copper  in  concentrations  up  to 
15  at.  %  nickel.   These  observations  provide  evidence  that  a  virtual - 
bound-state  model  applies  to  these  alloy.'.-,  in  agreement  with  results 
obtained  by  photoemission  studies. 
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When  zinc  is  alloyed  to  copper,  the  fundamental  edge  nt  2.1  eV 
of  copper  was  shifted  to  higher  energies;  on  the  other  hand,  the  struc- 
ture at  about  4.0  and  4.4  eV  was  shifted  to  lower  energies.   This  behav- 
ior was  interpreted  using  existing  band  structure  calculations  for 
alpha-copper  zinc  alloys.   The  agreement  between  theory  and  experiment 
was  satisfactory. 

Results  for  alpha-copper-aluminum  alloys  were  similar  to  those 
for  the  copper-zinc  system. 
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CHAPTER   I 


INTRODUCTION 


In  optical  studies  of  metals  and  alloys  the  material  is  probed 
with  electromagnetic  radiation  in  the  spectral  region  from  the  infrared 
to  tiie  ultiinviolet.   Light  reflected  from  a  metal  surface  results  from 
interaction  of  electrons  in  the  solid  with  the  incident  radiation. 
The  reflected  radiation  provides  information  on  the  energy  distribu- 
tion of  electrons  within  the  bands,  and  particularly  on  the  optical 
excitations  involved  in  these  interactions. 

Knowledge  of  the  band  structure  of  metals  has  improved  consider- 
ably in  recent  years.   The  energy  distribution  of  electrons  within  the 
bands  i.s  in  general  understood  and  in  many  instances  con  be  calculated. 
Measurements  leading  to  the  knowledge  of  the  band  structure  were  at 
first- possible  only  for  electrons  near  the  Fermi  surface.   However  with 
optical  techniques,  which  excite  optical  transitions  over  a  region 
covering  a  few  electron  volts  above  and  below  the  Fermi  surface,  the 
possibilities  have  been  greatly  extended. 

The  increasing  awareness  of  the  importance  of  understanding 
the  electron  structure  of  solids  in  metals  research  can  perhaps  be 
emphasized  by  quoting  the  words  of  J.  J.  Oilman  and  W.  A.  Tiller  [1]. 
The  progress  made  in  the  knowledge  of  atomic  and  electronic  structure 
of  metals 


will  have  an  important  Influence  on  understanding  and  predict- 
ing metallurgical  reactions  and  processes,  and  Is  now  In  a  form 
that  allows  It  to  be  Incorporated  Into  the  general  education  of 
metallurgists  ...  as  critical  materials  become  more  special- 
ized, their  Invention  and  characterization  will  depend  as  much 
on  the  control  of  their  electronic  structure  as  on  the  tradi- 
tional control  of  chemical  composition  and  microstructure  [1,  p.  v] , 


Optical  properties  of  solids  received  considerable  attention 
during  the  1950' s.   This  was  in  part  due  to  the  advances  in  materials 
preparation  technology  and  to  theoretical  work  leading  to  detailed 
knowledge  of  the  band  structure  of  solids. 

Conventional  techniques  to  measure  the  optical  constants  of 
solids  use  the  method  known  as  elllpsometry.   For  an  absorbing  medium, 
such  ns  metol ,  the  determination  of  the  optical  constants  Involves  the 
measurements  of  two  parameters  for  each  frequency  for  a  given  incident 
radiation.   If  lineorly  polarized  light  impinges  obliquely  on  a  metal, 
the  reflected  light  will  in  general  be  elllptically  polarized.   It  is 
possible  to  determine  the  optical  constants  at  a  given  wavelength  by 
measuring  the  state  of  polarization  of  the  reflected  light,  that  is, 
the  orientation  and  axes  ratio  of  the  polarization  ellipse.   A  number 
of  techniques  have  been  devised  to  perform  this  type  of  measurements 
but  they  are  restricted  to  the  spectral  region  for  which  polarizing 
equipment  and  wave  plates  are  available.   Moreover,  they  are  limited 
in  resolution  and  the  measurements  are  often  tedious  and  time  consum- 
ing, particularly  when  the  optical  constants  are  to  be  obtained  over 
a  range  of  wavelengths. 

More  recently  optical  studies  have  been  carried  out  at  normal 
or  near  normal  Incidence,  in  conjunction  with  relations  by   which  it  is 


possible  to  calculate  one  optical  constant  if  the  other  is  measured 
over  a  sufficiently  wide  spectral  range.   The  common  practice  is  to 
measure  the  reflectivity  and  to  determine  the  optical  constants  (usu- 
ally the  real  and  imaginary  parts  of  the  dielectric  constant)  from 
these  reflectivity  data.   With  normal  incidence  reflection  techniques, 
studies  were  extended  in  the  late  1950' s  and  early  1960's  to  photon 
energies  covering  the  infrared,  visible  and  ultraviolet.   By  the  mid 
I960' s  the  electronic  stx'ucture  and  the  optical  constants  of  a  large 
number  of  materials  had  been  investigated  and  tentatively  Interpreted. 
From  these  istudies  the  existence  of  the  so-called  critical  points  asso- 
ciated with  relatively  narrow  regions  in  frequency  were  discovered,  which 
produce  pronounced  structure  in  the  optical  constants.   ITie  Identifica- 
tion of  this  structure  with  particular  optical  transitions  today  constl- 
tutes  part  of  the  main  effort  in  the  area  of  optical  properties  of 
solids.    j 

Optical  properties  regained  Interest  in  ?.965  when  Seraphin  and 
Hess  [2J  introduced  o  modulation  technique  which  prompted  a  series  of 
new  experiments. 

The  optical  structure  associated  with  critical  points  can  be 
greatly  enhanced  by  means  of  modulation  or  derivative  techniques. 


The  general 


idea  behind  this  method  is  as  follows.   The  dielectric 


constants  near  a  three-dimensional  critical  point  can  be  expressed  as 

i 
e({u)  =  b(uj -uu  )«  +  constant,  (a) 

e 

where  cu  represents  the  critical  frequency,  m   is  the  frequency  of  the 
incident  radiation,  and  b  is  a  parameter.   The  constant  background 


In  Relation  (a)  may  be  much  larger  than  the  singular  part,  and  there- 
fore the  singularity  may  be  lost  in  the  noise  of  the  background.   It 
is  therefore  convenient  to  measure,  instead  of  e  (or  some  other  optical 
parameter) ,  the  derivative  of  e  with  respect  to  some  parameter  p. 
The  constant  background  is  thus  eliminated.   The  derivative  becomes 


de   b,     .-i   d.     -    db.     ,4  , ^ . 

_  =  -to  -u;^)  5  .  -to  -u,^)  +  -to  -C0g)1^  .  (b) 


It  can  be  seen  that  the  first  term  in  Relation  (b)  exhibits  a   singular- 
ity at  the  critical  frequency  u)   and  may  therefore  be  easily  detected. 
The  second  term  containing  to  -u)  )^  on  the  other  hand  can  be  neglected 

near  w  . 
g 

Several  parameters  p  can  be  used  to  modulate  the  band  struc- 
ture.  Seraphin  and  Hoss  [2]  imposed  an  alternating  electric  field  on 
a  semiconductor  to  induce  changes  in  the  reflectivity  (electroref lect- 
ance) .   Engeler  etal.  [3]  varied  an  external  uniaxial  stress  in  var- 
ious crystallographic  directions  during  reflectivity  measurements 
(piezoreflectance) .   Berglund  [4]  and  Hanus,  Feinleib  and  Scouler  [5] 
modulated  the  temperature  of  a  specimen  on  optical  absorption  measure- 
ments.  It  is  also  possible  to  obtain  the  first  derivative  of  the 
spectral  reflectivity  by  imposing  a  small  periodic  variation  on  the 
frequency  of  the  incident  radiation  (frequency  modulation).   Each  of 
these  modulation  methods  has  its  specific  advantages  and  disadvantages. 

Electroreflectance  has  been  used  successfully  in  the  study  of 
the  band  structure  of  semiconductors;  but  it  does  not  find  application 
In  metals  because  of  the  difficulty  of  maintaining  an  electric  field 


within  a  conductor.   Temperature  modulation  is  used  to  study  the  influ- 
ence of  phonons  on  interband  transitions.   Modulation  of  uniaxial 
stress  upon  a  sample  changes  the  symmetry  of  the  crystal  and  this  is 
found  useful  in  determining  the  symmetry  of  the  critical  points. 
Frequency  modulation  enhances  the  structure  in  the  reflectivity  but  It 
has  the  disadvantage  that  it  also  reveals  any  structure  present  in  the 
light  source. 

The  investigations  herein  reported  use  another  type  of  modula- 
tion:  the  variation  of  the  composition  of  alloys.   The  addition  of  an 
alloying  element  to  a  metal  as  a  rule  changes  its  band  structure. 
Compositional  modulation  provides  therefore  a  means  to  study  the  evolu- 
tion of  the  band  structure  of  the  alloy  as  its  oomposition  is  changed. 

The  differential  reflectometer  which  was  developed  in  this  work 
has  proved  to  be  potentially  a  powerful  tool  for  the  investigation  of 
the  optical  properties  of  solids  and  the  study  of  surface  reactions. 
Generally  speaking  the  instrument  is  capable  of  detecting  minute  reflec- 
tivity  differences  (a  fraction  of  1  percent  at  any  given  wavelength) , 
The  change  in  reflectivity  might  be  due  to  alloying,  a  surface  reaction, 
a  difference  in  doping  in  a  semiconductor,  a  change  from  crystalline  to 
amorphous  structure,  a  transformation  (short  range  order,  allotropic, 
etc, ) ,  or  any  phenomenon  which  brings  about  a  change  in  the  electronic 
structure  of  the  solid  with  which  a  change  in  optical  parameters  is 
associated. 

Because  of  computer  elements  incorporated  into  the  instrument, 
the  data  are  processed  instantaneously.   This  permits  automatic 


recording  of  a  normalized  change  in  reflectivity  AR/R  as  a  function 
of  wavelength.   Thlis  normalization  is  extremely  useful,  for  it  makes 
the  recording  of  the  data  independent  of  source  instability,  spectral 
sensitivity  of  detector,  and  intensity  of  light  source. 

Several  models  have  been  suggested  to  explain  the  effect  of 
alloying  on  the  electronic  structure.   Notably  among  these  are  the 
virtual-bound  state  and  rigid-band  models.   In  general  different  models 
lead  to  different  predictions  about  the  effects  of  alloying  on  the  band 
structure.   Therefore  a  study  of  the  changes  in  reflectivity  upon  adding 
solute  atoms  to  a  metal  can  provide  a  critical  test  for  the  validity  of 
current  models. 

Chapter  II  presents  an  outline  of  the  theory  from  which  inter- 
pretation of  the  optical  properties  of  metals  is  made.   The  optical 
constants  are  defined,  and  their  relation  to  electrical  properties  of 
the  solid  is  established.   A  brief  discussion  of  early  models  and  their 
relevance  to  modern  developments  is  also  given.   The  chapter  proceeds 
with  a  description  of  the  optical  excitation  processes  involved  in  the 
interaction  of  light  in  solids.   The  remainder  of  that  chapter  is 
devoted  to  a  review  of  some  of  the  more  relevant  achievements  in  the 
area  of  optical  properties  of  metals. 

The  description  of  a  new  apparatus  developed  for  optical 
measurements  by  compositional  modulation,  including  the  circuit  used 
to  process  the  data,  is  given  in  Chapter  III. 

In  Chapter  IV  the  results  obtained  by  this  technique  on  copper 
based  alloys,  with  additions  of  nickel,  zinc  and  aluminum,  are  pre- 
sented and  discussed.   The  main  effort  is  directed  towards  the 


correlation  of  the  experimental  findings  with  other  studies  and  to 
current  theoretical  models  of  alloying.   In  the  case  of  copper-nickel 
alloys  it  is  found  quite  conclusively  that  a  virtual-bound  model  applies 
to  copper-rich  alloys,  in  agreement  with  photoemission  studies  of  these 
alloys.   Measurements  on  the  copper-zinc  system  are  compared  with  band 
structure  calculations  on  these  alloys.   The  correlations  are  found  to 
be  satisfactory  in  some  respects  and  not  so  satisfactory  in  others. 
Results  obtained  on  copper-aluminum  alloys  are  discussed  in  terms  of 
the  evolution  of  its  structure  as  a  function  of  composition.   It  is 
suggested  that  this  system  probably  conforms  to  a  rigid-band  model. 


CHAPTER  II 
OPTICAL  PROPERTIES  OF  METALS 

2. 1  Introduction 

From  the  point  of  view  of  the  propagation  of  electromagnetic 
waves  in  an  absorbing  medium,  the  optical  properties  of  the  medium  are 
given  in  terms  of  two  parameters:   the  index  of  refraction,  n,  and  the 
absorption  coefficient,  k.   The  two  •  arameters  may  be  defined  by  the 
following  equations.   First, 

n  H  £  ,  (1) 

where  c  is  the  velocity  of  light  in  vacuum  and  v  is  the  phase  velocity 
in  the  medium.   This  is  the  definition  usually  given  for  the  refractive 
index  of  a  dielectric  substance.   'Hie  absorption  coefficient  k  is 
defined  by  the  rate  of  decrease  of  amplitude  E  of  the  electromagnetic 
wave  as  it  progresses  through  the  medium 


2tt 
E  =  E   exp  (-  -^  k  2)  .  (2) 

o         \ 


Here  z  is  the  coordinate  in  the  direction  of  propagation,  E   is  the 
amplitude  at  z  =  0,  and  \   is  the  wavelength  of  the  light  in  vacuum. 
The  name  "absorption  coefficient"  for  k  is  perhaps  not  the  best  choice, 
since  it  is  not  related  to  the  absorption  of  radiation  inside  matter  in 
a  very  direct  or  intuitive  way.   Extinction  coefficient  would  be  a 

8 


better  name  but  unfortunately  this  designation  Is  reserved  for  a  dif- 
ferent parameter  defined  as  k/n.   The  term  absorption  coefficient  will 
be  used  here  since  It  has  become  the  standard  custom. 

Together  these  two  parameters  form  the  complex  Index  n  of  the 
medium.   This  is  defined  by  the  equation 

n  =  n  -  ik  .  (3) 

The  term  "complex  index"  for  the  quantity  defined  by 
Equation  (3)  is  an  approprl  te  one.   If  the  Index  in  the  equation 

E  =  E  exp  I  iu)(t  -^  z) J  (4) 

for  transpai^ent  media,  is  replaced  by  the  complex  Index  n,  the  result- 
ing equation, 


E  =  E^  exp  (-  ^  k  /.  )  oxp  [  IcuCt  "  J  z)J  , 


(5) 


For  example, 


c  optics  to  be  converted  easily  to  the  metallic  case. 

2 


has  the  properties  required  by  the  defining  Equations  (1)  and  (2). 
This  concept  is  convenient  because  it  enables  some  of  the  formulas 
of  dielectric 

the  reflectivity  of  a  dielectric   is    | (n  -  l)/(n + 1)  |      . 
It  can  be  shown   [6]    that   the  reflectivity  U  of   a  metal   is  found  by 
inserting  the  complex  index  into  the  corresponding  expression  for   the 
reflectivity  of   a  dielectric,    namely, 

2  2 

(n-1)      +  k  ,    . 

R  =  J- g-     .  (6) 

(n+  1)      +  k 
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The  conversion  of  optical  formulas  to  the  metallic  case  In  this 
way  provides  little  insight  into  the  physics  involved  and  so  gives  no 
aid  in  understanding  or  interpreting  the  equations.   It  is  preferable, 
therefore,  to  treat  the  absorbing  case  directly.   Transparent  media  can 
then  be  considered  as  a  special  case  In  which  k  =  0. 

2.2  Relation  Between  the  Optical 
and  Electric  Parameters 

The  classical  theory  of  the  optical  properties  of  solids  is 

based  upon  Maxwell's  electromagnetic  equations.   When  the  solid  Is  n 

conductor  it  is  cliaracterized  by  a  conductivity  a  and  a  ditlectric 

constant  e.   The  permeability  of  the  medium  is  assumed  unity.   For  an 

uncharged  medium  the  equations  take  the  form 


div  MI  =  0      curl  H  =  ^  E  +  i  ^ 
*^  c      c  ot 

(7) 
div  c  E  =  0      curl  E  =  -  -  ^ 


where  the  Gausdlan  system  of  units  is  used.   The  equations  as  written 
are  supplemented  by  Ohm's  law,  J  =  a  E. 

For  a   homogeneous,  isotropic  (or  cubic)  medium  the  dielectric 
constant  and  the  conductivity  are  scalars  Instead  of  tensors.   The 
magnetic  field  H  can  be  eliminated  between  the  two  curl  equations  and, 
using  the  fact  that  div  E  =  0,  the  wave  equation  can  be  derived,  namely, 

-•      2-* 

-;   4TTa  3e   e  S  E  ,Q. 

^  =  =  —  ^  +  -2  rr  •  ^»^ 

c       c  ot 
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This  is  the  differential  equation  that  governs  the  propagation  of 
electromagnetic  waves  in  a  conducting  medium.   It  can  be  seen  that  the 
term  containing  the  first  derivative  of  E  is  a  dissipative  term.   It  can 
be  shown  [7]  that  In  a  steady  state  situation  no  longitudinal  electric 
or  magnetic  field  (other  than  uniform  fields)  can  exist  in  a  conductor 
which  satisfies  Equation  (8).   It  is  possible  therefore  to  write 
Equation  (8)  in  a  one-dimensional  form.   With  the  direction  of  vibra- 
tion of  the  electric  field  along  the  x-axis  and  the  direction  of  propa- 
gation  aloiig  the  z-axis,  the  equation  takes  the  form 


B^E    .   Se     „  b\ 

X    4TTCT    X     €      X 

oz     c        c   ot 

A  plane-wave  solution  of  (9)  of  angular  frequency  U)  will  then  have 
the  form 

E^  =  E^  exp  |^iuj(t  -B  z) J  ,  (10) 


where  the  index  of  refraction  is  complex.   Substitution  of  solution  (10) 
into  the  wave  equation  (9)  yields 


-2       4-na 
n  =  e  -  i  


(n-ik)2=  ^-i^.  (11) 


which  establishes  a  relation  between  the  optical  constants  (n,k)  with 
the  electric  properties  of  the  solid  ie,a)   for  the  particular  frequency tu. 
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It  will  be  observed  from  Equation  (11)  that  it  is  possible  to 
introduce  a  complex  dielectric  constant 


c  =  e  -  i  ins  .  (12) 


Here  the  real  part  of  e  is  the  conventional  (actual)  dielectric  constant, 
while  the  imaginary  part  Includes  the  disslpative  effects  associated  with 
the  conductivity.   This  is  a  convenient  generalization  of  the  dielectric 
constant,  since  it  allows  the  relation 

r?   =  e  (13) 

for  transparent  substances  to  be  extended  to  the  metallic  case  with 
n  and  e  substituted  by  the  complex  quantities  introduced  above. 

It  appears  convenient  to  look  at  the  complex  nature  of  the 
dielectric  constant  by  considering  the  curl  H  equation  in  (7).  It  can 
be  seen  that  the  contribution  to  the  magnetic  field  comes  from  conduc- 
tion currents  through  the  conductivity  and  also  from  displacement  cur- 
rents  arising  from  the  time  derivative  of  the  displacement  vector  D. 
For  a  particular  angular  frequency  (u  the  time  dependency  of  the  fields 
will  be  of  the  form 

E  =  E   exp  (itut)  , 
o 

and  the  time  derivative  will  be 


s=-? 
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With  this  last  expression  it  is  therefore  possible  to  formally  trans- 
form the  right-hand  side  of  the  curl  H  equation  into  terms  containing 
only  the  time  derivative  of  the  electric  field,  that  is, 

curl  H  =  -  (^c  -  i  _j  ^  .  (14) 

The  two  terms  in  parentheses  represent  the  same  complex  dielectric 
constant  introduced  in  Equation  (11). 

In  a  similar  way  it  is  also  possible  to  generalize  the  conduc- 
tivity by  converting  the  right-hand  side  of  Equation  (14)  into  terms 
containing  only  the  field  E.   The  resultant  relation  is 

c„rlH=f  (<,*1^)e,  a6) 

where  the  generalized  conductivity  can  be  defined  as 


a  =  a  +  i  ^  .  (16) 


It  should  be  emphasized  that  this  is  just  a  convenient  general- 
Ization  of  the  parameters  involved  in  the  Interaction  of  electromagnetic 
waves  with  matter.   Whichever  formalism  is  used,  the  relation  between 
the  optical  and  electrical  parameters  is  the  one  given  by  Equation  (11), 
The  aim  of  a  theory  of  the  optical  properties  is  to  obtain  relations 
for  n  and  ki  or  equivalently  for  e  and  a  based  on  a  microscopic  model 
of  the  conductor.   If  a  particular  model  of  the  conductor  leads  to  an 
equation  relating  the  dielectric  constant  to  fundamental  parameters 
such  as  the  mass  and  charge  of  the  electron  and  if  the  dielectric 
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constant  Is  a  complex  number,  then  the  imaginary  part  can  be  associated 
with  conduction  currents,  that  is,  with  the  conductivity.  On  the  other 
hand,  if  one  is  dealing  with  a  complex  conductivity,  it  can  be  asserted 
that  its  imaginary  part  is  taking  care  of  the  polarization  current. 

2.3  Optical  Excitation  Processes 

The  distinction  between  metals  and  insulators  is  made  on  the 
basis  of  dc  electricnl  properties.   Upon  the  application  of  a  steady 
electric  field,  a  metal  conducts,  whereas  an  Insulator  does  not.   This 
distinction  becomes  increasingly  vague  when  the  applied  field  is  fre- 
quency dependent,  particularly  when  the  frequency  falls  into  the  optical 
range,  that  is,  into  the  visible  and  ultraviolet.   In  this  range  the 
electrical,  or,  more  appropriately,  the  optical  properties  of  metals  and 
insulators  have  much  in  common.   It  seems  appropriate  therefore  to 
start  this  section  with  a  discussion  of  some  elementary  properties  which 
apply  to  both  metals  and  dielectrics. 

The  earliest  useful  model  for  the  electrical  and  optical  prop- 
erties of  metals  was  proposed  by  Drude  early  this  century.   In  this 
model  it  is  supposed  that  the  electrons  form  a  gas  of  free  particles 
which  can  respond  to  the  application  of  external  fields  by  producing 
currents.   The  ionic  cores  present  in  the  solid  are  completely  neglected. 
The  only  statistics  known  at  that  time  was  that  of  Boltzmann.   It  was 
later  discovered  that  through  the  substitution  of  Fermi  statistics  for 
Boltzmann  statistics,  Drude' s  picture  could  be  transformed  into  a  very 
realistic  and  useful  model  for  a  metal.   Even  in  its  original  form, 
Drude' s  model  is  already  adequate  for  the  discussion  of  many  aspects 
of  the  optical  properties  of  metals. 
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Another  simple  picture  of  an  Insulator  was  proposed  by  Lorentz. 
Here  the  electrons  are  envisioned  as  tied  to  certain  equilibrium  posi- 
tions by  electrical  forces  arising  from  the  ionic  cores.   Since  the 
electrons  are  not  free  to  wander  around,  the  crystal  can  have  no  dc 
conductivity.   On  the  other  hand,  polarization  effects  in  the  presence 
of  an  nc  electric  field  are  quite  possible. 

In  the  light  of  sophisticated  modern  developments  in  solid 
state  physics,  it  is  perhaps  remarkable  that  such  simple  ideas  still 
have  releyance  in  the  interpretation  of  the  optical  properties  of 
solids.   If  one  can  understand  qualitatively  how  electrons  described 
by  these  models  respond  to  light,  or  equivalently  to  an  ac  electric 
field,  then  a  synthesis  and  generalization  of  these  ideas  will  come 
close  to  a  true  description  of  the  complex  nature  of  the  interaction 
of  light  wilth  the  electrons  in  a  real  solid. 

2.3.1  Classical  Oscillators 

The  response  of  the  electrons  in  a   Lorentz  Insulator  to 
a  uniform  oscillating  electric  field,  E(t)  =  E(u))  exp  (-icot)  ,  having 
a   single  Fourier  component  corresponding  to  frequency  u),  is  described 
by  the  equation  of  motion 


r 


2 
— I  ■^?dT  +  "^  ^)~  6E(u))  exp  (-iiut)  .  (17) 


dt 


2    2 
Here  d  x/dt  is  the  acceleration  term,  l/T(dx/dt)  is  a  damping  terra 

proportional  to  the  velocity  and  characterized  by  the  relaxation  time  T 

2      I 
and  cu  x  is  associated  with  the  restoring  force  of  the  oscillator  whose 

s 
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natural  frequency  is  uj  .   The  right-hand  side  of  the  equation  repre- 
sents the  driving  term  due  to  the  electric  field  with  e  and  m  being 
the  electric  charge  and  the  mass  of  the  electron,  respectively.   Since 
the  equation  is  linear  in  the  displacement  x,  it  can  be  assumed  that 
x(t),  and  all  electromagnetic  field  quantities,  have  the  same  time 
dependence  as  E(t);  that  is,  x(t)  =  x(u))  exp  (-iujt).   The  substitution 
of  x(t)  transforms  Equation  (17)  into  a  linear  algebraic  equation  for 
x(u))  that  can  be  solved  readily. 

The  vibrations  of  all  the  electrons  in  a  unit  volume  of  the 
solid  result  in  an  induced  polarization  given  by 

P(t)  =  P((u)  exp  (-icut)  (18) 

whose  amplitude,  P(u))  =enx(u)),  is  proportional  to  the  displacement 
x((ju)  and  the  concentration  n  of  electrons.   Using  the  value  x((ju)  that 
comes  from  the  solution  of  Equation  (17),  P(u))  is  found  to  be 


2 

T,/  N       "e  E(u))/m  ,,_v 

P(u))  = g 2 •  ^^^' 

U)  -10  +  iu)/T 

s 


The  frequency-dependent  dielectric  constant  describing  the  response 
of  the  system  to  an  external  field  E(iu)  exp  (-iiwt)  having  frequency  to 
is  related  to  the  electric  field  and  the  polarization  by  the  expression 

e(io)  E(t'j)  =  E(u))  +  4ttP(u))  (20) 

using  cgs  units.   Combininj  (19)  and  (20)  the  complex  dielectric 
constant  Is  obtained 
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U) 

e(aj)  =  1  -  — — (21) 

CO  -U)  +  iU)/T 

s 

2    8 

where  (u  =  (4Tme  /m)   is  the  so-called  plasma  frequency  associated 
with  electrons. 

It  should  be  emphasized  that  the  dielectric  constant  is  here 
a  complex  quantity.   The  real  part  is  e  as  It  appears  in  Maxwell  equa- 
tions and  the  imaginary  part  includes  the  dissipative  effects  asso- 
ciated with  the  conductivity.   In  fact,  the  model  lending  to  Equa- 
tion (21)  contains  both  the  Drude  model  and  the  Lorentz  oscillntors. 
Equation  (21)  does  not  contain  the  conductivity  explicitly,  but  if  the 
imaginary  part  is  separated  out,  the  conductivity  can  be  shown  readily 
to  be 

CT  =  a ^  (22) 

°  1  +  (mr 

2 
where  a  =  ne  T/m  is  the  dc  conductivity.   It  can  be  observed  that. 

o  ''  ' 

by  introducing  a  relaxation  time  T  ^  0  into  the  equation  of  motion  (17), 
a  dissipative  contribution  comes  in,  that  is,  the  conductivity  of  the 
solid.   Furthermore,  the  real  and  imaginary  part  of  the  dielectric 
constant  (?1)  can  be  related  to  the  index  of  refraction  n  and  the 
absorption  coefficient  k  through  the  formula  (11)  derived  in  Section  2.2. 
Thus,  the  complex  dielectric  constant  and  the  parameters  n  and  k  contain 
the  same  physical  information.   For  purposes  of  theoretical  descrip- 
tions, howeyer,  the  dielectric  constant  formulation  is  more  convenient 
because  it  can  be  related  in  a  more  direct  way  to  a  microscopic  model 
of  the  solid. 
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The  behavior  of  the  function  described  in  Equation  (21)  iu  shown 

in  Figure  2.1  where  e,  and  e  represent  the  real  and  Imaginary  parts, 
1      A 

respectively.   For  an  insulator,  in  which  the  restoring  force  ia  finite, 
the  dielectric  constant  is  that  of  an  oscillator,  Flgiire  ?..1A.   On  the 
other  hand,  for  a  metal,  the  Drude  model  can  be  recovered  simply  by 
equating  the  natural  frequency  oo  to  zero,  muanliig  that  the  electrons 
are  free  to  move  through  the  crystal.   The  characteristic  behavior  for 
the  real  and  imaginary  parts  of  the  dielectric  constant  is  now  shown  in 
Figure  2. IB;  e  and  c  become  positi'ely  and  negatively  large,  respec- 
tively, with  diminishing  frequency.   Thus,  it  is  possible  to  character- 
ize the  electrons  in  a  Drude  metal  ns  oscillators  of  zero  natural 
frequency. 

2.3.2  Energy  Bands  in  Solids 

As  mentioned  earlier,  Drude' s  model  of  a  metal  completely 
ignores  the  presence  of  the  ionic  cores.   On  the  other  hand,  the  Lorentz 
model  of  an  insulator  takes  another  extreme  position:  the  electrons  are 
bound  to  the  positive  ions  with  elastic  forces  and  can  only  oscillate 
about  their  equilibrium  positions  but  cannot  move  around  inside  the 
crystal.   The  situation  in  a  real  metal  has  features  of  both  the  free- 
electron  model  and  the  Lorentz  insulator.   Some  electrons  are  tightly 
bound  to  the  positive  nuclei  and  are  little  affected  by  the  fact  that 
the  atoms  are  part  of  a  solid,  while  other  electrons  move  almost  freely 
inside  the  crystal. 

The  model  that  emerges  when  the  interaction  between  the  elec- 
trons and  the  periodic  potential  of  the  crystal  is  taken  into  account 
is  considerably  more  complex  than  the  simple  models  described  above 
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Figure  2. 1  Behavior  of  the  function  described  in  Equation  (21)  where 
€-,  and  62  represent  the  real  and  Imaginary  parts,  respec- 
tively,   A  -  Insulator  (U)  ^  0).   B  -  Metal  (U)  =  0). 
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and  Its  development  requires  sophisticated  techniques  of  quantum  mechan- 
ics.  The  description  of  the  details  leading  to  the  band  theory  of 
solids  Is  beyond  the  scope  of  this  dissertation.   However,  It  seems 
appropriate  to  present  here  some  of  the  underlying  principles  of  the 
theory  as  well  as  some  of  Its  prominent  features. 

The  electron  In  an  Isolated  atom  can  only  have  certain  pre- 
scribed energy  states  and  as  a  result  the  electron  can  change  Its 
energy  only  by  discrete  amounts.   When  the  atoms  are  put  together  to 
form  a  crystalline  solid,  the  levels  of  the  Individual  atoms  spread  out 
Into  bands.   This  phenomenon  can  be  regarded  as  a  consequence  of  the 
Paull  exclusion  principle  which  comes  into  effect  when  the  wave  func- 
tions of  the  various  atoms  overlap  In  the  solid.   The  amount  of  spread 
Is  related  to  the  extent  of  overlapping  of  the  wave  functions.   The 
energy  states  Inside  the  crystal  can  thus  be  envisioned  as  a  series  of 
bands  with  a  quaslcontlnuous  distribution  of  states  within  each  band: 
narrow  bands  for  the  more  tightly  bound  electrons  with  the  width  becom- 
ing progressively  larger  for  higher  energies  and  eventually  overlapping 
with  each  other  at  sufficiently  high  energies.   The  energy  difference 
between  electronic  states  of  an  Isolated  atom  become  the  "energy  gaps" 
in  a  crystal.   There  are,  therefore,  energy  regions  not  allowed  or 
forbidden  gaps. 

In  a  solid,  energy  levels  are  generally  described  in  terms  of 
a  wave  vector  k  (|k|=  2tt/X)  and  a  band  Index.   In  band  calculations 
one  of  the  tasks  of  quantum  mechanics  is  to  find  a  relation  between  the 
energy  E  and  the  vector  k  for  each  band.   For  free  electrons  this  is 
a  simple  relation,  namely. 
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2 
E  =  2_  k^  (23) 

2m 

where  h  =  h/2TT  (h  =  Planck's  constant)  and  m  is  the  electronic  mass. 
When  the  electron  moves  in  the  periodic  potential  of  the  lattice,  how- 
ever, Equation  (23)  is  not  valid.   In  dealing  with  the  motion  of  elec- 
trons in  a  crystalline  solid  it  is  convenient  to  work  with  the  so-called 
reciprocal  lattice  or  reciprocal  space.   A  consequence  of  the  periodic- 
ity of  the  lattice  in  ordinary  space  is  that  the  reciprocal  space  is 
also  periodic.   When  the  wave  vector,  which  is  o  vector  in  reciprocal 
space,  approaches  one  of  the  vectors  of  the  reciprocal  lattice,  the 
condition  for  Bragg  reflection  is  satisfied.   Bragg  reflection  occurs 
also  for  electron  waves  in  crystals  and  leads  to  energy  gaps.   When 
the  Bragg  condition  is  satisfied  for  electrons  of  wave  vector  k,  waves 
travelling  in  one  direction  are  soon  Bragg -reflected  and  then  travel 
in  the  opposite  direction,  leading  to  standing  waves.   In  other  words, 
no  travelling  waves  of  the  wave  equation  can  exist  and  therefore  devia- 
tions from  the  simple  parabolic  behavior  of  Equation  (23)  are  expected. 
As  an  example  for  a  real  metal,  the  bands  along  certain  symmetry  direc- 
tions are  shown  In  Figure  2.2  for  aluminum.   As  it  happens,  aluminum 
has  a  band  structure  which  deviates  very  little  from  the  free-electron 
bands  shown  in  the  figure  by  the  dashed  curves.   The  free-electron 
bands  are  simply  obtained  by  reducing  the  interaction  of  the  electrons 

with  the  periodic  potential  to  a  negligible  amount  while  keeping  the 

i 
symmetry  properties  of  the  lattice.   It  is  customary  to  present  the 

bands  along  certain  directions  of  particular  symmetry  in  k  space.   These 

symmetry  directions  and  points  are  shown  in  Figure  2.3.   Also  shown  in 
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Figure  2.2   Energy  bands  along  two  symmetry  directions  for  aluminum. 
Solid  curves  represent  calculated  energies;  dashed  curves 
show  the  free-electron  energy  bands.   a  -  Intraband  tran- 
sition,  b  -  Interband  transition. 


Figure  2.3   Brillouin  zone  for  the  face-centered  cubic  structure  show- 
ing various  points  and  lines  with  symmetry  properties. 
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this  figure  is  the  fundamental  unit  cell  of  the  face-centered  cubic 
reciprocal  lattice  also  called  the  Brillouin  zone.   One  practical 
importance  of  this  cell  is  that  the  wave  vectors  corresponding  to 
independent  solutions  of  the  wave  equation  are  all  contained  in  this 
unit  cell.   That  is,  the  wave  equation  needs  only  be  solved  for  k 
vectors  contained  within  the  Brillouin  zone.   Although  Figure  2.2 
represents  the  energy  levels  for  a  single  electron  only,  the  electrons 
in  a  solid  can  be  regarded,  to  a  good  approximation,  as  independent  and 
the  diagram  can  be  applied  to  such  systems.   In  this  independent  par- 
ticle model  the  electrons  know  of  one  another' s  existence  only  througli 
the  exclusion  principle,  which  provides  that  no  more  than  two  electrons, 
one  of  each  spin,  can  occupy  a  given  state.   This  means  that  at  low 
temperatures  all  levels  below  the  Fermi  level,  indicated  in  Figure  2.2 
by  the  horizontal  line,  are  filled  and  all  others  are  empty. 

2.3.3  Excitation  Processes 

in  this  section  the  three  relevant  optical  excitation 
processes  in  metals  will  be  described.   Intraband  transitions  are 
analogous  to  the  optical  excitations  in  a  free-electron  metal.   As  indi- 
cated in  Figure  2.2,  intraband  transitions  involve  the  excitation  of 
an  electron  by  a  photon  within  the  same  band.   Because  of  the  exclu- 
sion principle,  the  electron  must  go  from  an  initial  filled  state  to 
a  final  empty  state.   These  transitions  therefore  can  only  occur  in 
metals,  but  not  in  insulators,  in  which  all  bands  are  either  completely 
full  or  completely  empty.   By  contrast,  the  second  kind  of  excitation, 
the  interband  transition,  which  involves  two  bands,  is  common  to  both 
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insulators  and  metals.   This  process  has  its  analog  in  the  optical 
excitation  of  a  Lorentz  insulator  and,  in  fact,  reduces  almost  precisely 
to  that  for  the  very  artificial  case  of  a  solid  consisting  of  a  widely 
separated  array  of  two-quantum-level  atoms  whose  wave  functions  do  not 
overlap  appreciably.   In  the  absence  of  any  other  interactions  involv- 
ing the  electrons,  the  interband  transitions  can  be  shown  to  be  essen- 
tially vertical  in  k  space,  as  a  result  of  the  fact  that  the  total 
crystal  momentum  of  electron  plus  photon  must  be  conserved  during  the 
transition.   For  wavelengths  in  the  visible  and  ultraviolet  regions, 
the  momentum  of  the  light  wave  is  very  much  smaller  than  the  diameter 
of  a  Brillouin  zone,  and  therefore  the  wave  vector  of  the  electron  is 
essentially  unchanged  in  the  transition. 

Comparison  of  the  results  of  a  quantum  mechanical  derivation 
of  the  dielectric  constant  of  a  solid  with  the  elementary  results  out- 
lined above  shows  that  it  is  legitimate  to  regard  the  dielectric  con- 
stant as  a  sum  of  contributions  from  both  intraband  and  interband 
transitions.   Tliese  transitions  can  be  described  as  in  the  Drude 
free-electron  metal  plus  a  contribution  of  oscillators  at  each  k  con- 
necting vertically  filled  and  unfilled  states. 

Since  interband  transitions  are  possible  everywhere  In  the 
Brillouin  zone  and  for  all  energies  beyond  a  certain  threshold,  these 
transitions  might  be  expected  to  present  relatively  structureless 
contributions  to  the  dielectric  constant.   The  first  experimental 
evidence  that  this  is  not  the  case  came  in  connection  with  semiconduc- 
tor i.tudies  by  J.  C.  Phillips  [8]  where  very  large  contributions  from 
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certain  small  regions  In  k  space  were  found  that  result  In  pronounced 
structure  in  the  optical  parameters.   This  was  very  important,  since 
it  led  to  a  new  branch  of  an  old  discipline,  which  can  be  called  "band- 
structure  spectroscopy."  These  regions  in  k  space  associated  with 
strong  interband  transitions  contain  the  so-called  "critical  points." 
They  are  characterized  by  very  rapidly  varying,  often  large  density  of 
states  for  optical  transitions.   For  the  most  part,  these  occur  at 
points  of  high  symmetry,  points  at  which  the  band  structure  of  a  given 
material  is  most  likely  to  be  known.   Thus,  given  a  rough  idea  of  the 
band  structure,  it  is  possible  to  use  optical  techniques  to  obtain  some 
important  band  gaps,  thereby  providing  a  basis  for  better  calculations. 

The  third  kind  of  elementary  excitation  in  a  solid  is  the  so- 
called  plasma  or  collective  oscillation,  which  results  from  the  high 
density  of  electrons  that  are  present  in  a  metal  and  the  fact  that  they 
can  act  cooperatively  due  to  the  Coulomb  interaction  among  them.   This 
cooperative'  motion  may  be  regarded  as  a  fluid-like  vibration  of  the 
electrons  as  a  whole  through  the  solid,  which  produces  fluctuations 
in  the  electron  density.   Under  certain  conditions,  plasma  oscillations 
can  represent  normal  modes  of  the  entire  system.   This  means  that,  in 
the  absence  of  dlsslpative  effects,  once  such  oscillations  are  excited, 
they  do  not  decay  in  time.   Since  they  maintain  themselves,  an  Internal 
electric  field  E(u)  )  is  expected  at  the  plasma  frequency  U)p,  due  to  the 
oscillations  in  the  electron  density,  without  the  presence  of  an  exter- 
nal field  D(u)p).   Since 

D(u)p)  =  e(u)p)  E(u)p)  =  0 
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and  E((Up)  ^  0,    the  condition  for  plasma  oscillations  is  therefore 


e(u)p)  =   ej^((jUp)  +  iegCoUp)  =  0  .  (24) 


Thus,  strictly  speaking,  plasma  oscillations  can  exist  as  normal  modes 
of  the  system  only  if  both  the  real  and  imaginary  parts  of  the  dielec- 
tric constant  vanish  at  cu  .   In  practice,  however,  the  plasma  frequency 
will  be  well  defined  if  c  (oUp)  <  1  when  s^dUp)  =  0-   Since  e  repre- 
sents the  damping  of  the  plasma  resonance,  the  condition  e((j}_)  <   1 
Implies  that  the  damping  should  be  small. 

ilie  excitation  of  plasma  oscillations  can  be  described  by 
consideration  of  a  uniform  displacement  of  the  electron  gas  in  a  thin 
metallic  slab.   The  electron  gas  is  moved  as  a  whole  with  respect  to 
the  positive  ion  background.   A  displacement  of  amplitude  x  creates  an 
electric  field  E  =  4rmex  which  acts  as  a  restoring  force.   Therefore 
the  equation  of  motion  of  a  unit  volume  of  the  electron  gas  is 

2 
d  X       „     ^22  ,„_. 

nm  — g-  =  -  neE  =  -  4TTn  e  x  .  (25) 

dt^  "" 

Equation  (25)  is  recognized  as  the  equation  of  motion  of  a  simple 
harmonic  oscillation  of  frequency 


1 


=  [^J 


(26) 


This  expression  for  the  plasma  frequency  is  identical  with  that  appear- 
ing in  Equation  (21)  which  arose  in  quite  a  different  connection. 


27 


2.4  The  Interpretation  of  the  Optical  Properties 

The  linear  response  of  a  conducting  nonmagnetic  medium  to  an 
electromagnetic  wave  of  angular  frequency  tu  can  be  characterized  by 
a  frequency-dependent  complex  dielectric  constant  or  equivalently  by 
a  complex  index  of  refraction.   In  Section  2.2,  it  has  been  shown  that 
through  Maxwell's  equations  a  relation  exists  between  these  two  complex 
quantitiesL   From  the  point  of  view  of  the  interpretation  of  the  optical 
parameters  in  terras  of  a  theoretical  model,  the  complex  dielectric 
constant  plays  a  more  important  role  because  it  can  be  related  directly 
to  optical  excitation  processes.   Experimentally,  however,  the  dielec- 
tric constant  cannot  be  measured  directly  at  optical  frequencies  and 
therefore  reliance  has  to  be  made  on  reflectivity  measurements.   These 
measurements  are  based  on  Fresnel  formulas  for  the  reflection  and 
refraction  at  the  boundary  of  a  conducting  medium.   It  is  also  possible 
to  introduce  a  complex  reflectance.   This  is  defined  as  the  ratio  of 
the  complex  amplitude  of  the  reflected  to  that  of  the  incident  electric 
field.   The  complex  reflectance  r  =  r  exp  (10)  can  be  shown  [9]  to  be 
related  to  the  complex  Index  of  refraction  by  the  expression 


n  -  ik  -  1      19  ,„  , 

irriFTT=^^     •  <27) 


The  reflectivity  R  introduced  in  Section  5.1  is  obtained  immediately 

2 

by  R  =  r  .   Measurements  of  either  (n,k)  or  (r.S)  permit  therefore  the 

calculation  of  (^i.^g)  and  the  complete  description  of  the  optical  prop- 
erties in  the  appropriate  frequency  range  can  be  obtained.   Usually  the 
real  and  imaginary  parts  of  the  complex  index  of  refraction  are  measured 
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independently  for  each  frequency  of  light.   This  procedure,  however, 
is  long  and  tedious,  since  it  involves  measurements  at  different 
angles  of  incidence  for  each  frequency. 

Fortunately,  there  are  relations  between  the  real  and  imaginary 
parts  of  any  of  the  complex  quantities  describing  the  optical  properties. 
These  are  very  general  equations,  known  as  KramersHCronig  [10]  rela- 
tions, that  allow  the  calculation  of  one  parameter  at  any  given  fre- 
quency when  the  other  is  known  over  a  wide  range  of  frequencies. 
This  relationship  is  imposed  by  the  causal  nature  of  the  rehPonse  to 
an  electromagnetic  field.   That  is,  no  response  to  an  applied  field 
can  appear  before  the  field  is  applied.   For  the  dielectric  constant, 
this  relationship  is 


^l^'"^  -^  =  nPj   -72—2-*'  (28) 

o 

and  its  Inverse: 


(jU  -U) 

o 


where  P  designates  the  Cauchy  principal  value  of  the  integral.   For 
the  complex  reflectance  the  expression  Is: 


e.o.'ip  fjl|iiL>  d^'  .  (30) 


O  U)    -U) 
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The  use  of  the  Integral  (30)  requires  the  knowledge  of  the  reflectiv- 
ity for  all  frequencies,  in  principle.   Since  in  practice  its  measure- 
ment is  not  possible  for  all  frequencies,  an  explicit  assumption  for 
extrapolation  to  high  frequencies  is  necessary.   The  results  of  the 
method  are  more  reliable  if  the  measurements  of  R  are  extended  to 
higher  frejiuencies.   A  great  experimental  advantage  of  the  usage  of 
relation  (30)  lies  in  the  fact  that  only  the  reflectivity  at  or  near 
normal  incidence  has  to  be  measured.   Therefore,  the  method  using 
reflectivity  measurements  at  normal  incidence  has  been  widely  used  in 
recent  years. 

As  was  pointed  out  in  Section  2.3,  the  Interpretation  of  the 
optical  properties  of  metals  can  be  made  on  the  basis  of  a  dielectric 
constant  which  includes  contributions  from  free  electrons  plus  a  term 
associated  with  bound  electrons.   These  two  contributions  can  be 
thought  of  as  generalizations  of  the  simple  models  introduced  in 
Section  2.3;  the  free-electron  term  having  the  form  of  Drude  equation, 
whereas  the  bound  part  is  the  quantum  mechanical  counterpart  of  the 
Lorentz  equation  for  classical  oscillators.   Consequently,  the 


dielectric 
but ions  of 


constant  for  a  metal  is  composed  addltlvely  of  the  contri- 

(f) 
free  electrons  e   (cu) ,  and  Interband  transitions  or 


boimd  electrons  e   (uj) : 


e(cu)  =  e^*^(uj)  +  e^^^ui))    . 


Both  are.  In  general,  complex. 


30 


The  free-electron  contribution  is 


2 

U) 

C^^\(U)  =  1  -  —-J^l—-^  .  (31) 

^        U)(u)  +  i/T  ; 
c 

2       2 
Here  uj   =  dire  n  /m  Is  the  square  of  the  conduction  band  plasma  fre- 
Pa        c  a 

quency,  n  is  the  conduction  electron  density, 

±=        J-^        f  d\  i  Vl   E  (k)  (32) 

m    .  3. 2   J      3  k  c 
a       4ti  n  n 
c 

is  the  average  inverse  effective  mass  of  the  conduction  electrons,  the 

so-cnlled  "optical"  mass  defined  by  Cohen  [11].   The  two  parameters 

u)„  and  the  relaxation  time  t  characterize  the  conduction  or  free 
Pa  c 

electrons  in  the  Drude  theory. 

An  expression  for  the  interbnnd  processes,  and  valid  within 
the  fromework  of  the  random-phase  approximation  is 


Here  e  is  the  electron  charge,  m  the  free-electron  mass,  f «(k)  is  the 
Fermi  distribution  function  for  state  k  in  band  i, 

is  the  oscillator  strength, 
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IP*         a  3 

^i'l^   =  7-  J  "^'  (k,x)  P^  u^  (k,x)  d  X 


Is  the  momentum  matrix  element  for  direction  p,  integrated  over  the  unit 

cell  of  volume  v  ,  and 
a 

hU)^^'  =  E^  (k)-  E^/(k) 


is  the  energy  difference  between  bands  i   and  Z     at  the  point  k. 
T j,.t    is  an  interband  relaxation  time. 

The  above  formalism  has  been  developed  largely  by  Ehrenreich 

- 

and  his  co-workers  [12-15]  and  has  been  applied  with  a  great  deal  of 
success  to  the  interpretation  of  the  optical  properties  of  the  noble 
metals  as  well  as  the  metals  aluminum  and  nickel.   The  rest  of  this 
section  is  devoted  to  a  review  of  this  work  with  some  emphasis  on  the 
results  for  copper. 

Aluminum.  Figure  2.4  shows  the  reflectivity  R  of  aluminum  as 
a  function  of  the  energy  of  the  incident  radiation.  It  is  character- 
ized by  its  very  simple  optical  properties  that  make  this  metal  a  good 

example  of  a  free-electron  metal.   The  reflectivity  for  an  ideal  Drude 

■ 

free-electron  metal  stays  at  100  percent  as  long  as  e  is  negative. 
When  e  passes  through  zero  at  the  plasma  frequency ,  the  reflectivity 
falls  off  \^ery  rapidly.   It  can  be  seen  that  Al  appears  to  conform 
qualitatively  with  this  description  except  for  the  small  dip  near  1  eV, 
which  is  due  to  interband  transitions  and  therefore  represents  a 
distinct  departure  from  this  model. 
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Figure  2.4   The  spectral  dependence  of  the  reflectivity  of  aluminum. 
(From  Reference  14, ) 
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The  results  of  a  Kramers-Kronig  analysis  of  the  reflectivity 
are  shown  in  Figure  2.5.   It  can  be  observed  that  the  real  and  imagi- 
nary parts  of  the  dielectric  constant  near  zero  frequency  have  the 
behavior  characteristic  of  a  Drude  free-electron  metal.   Around  1  eV 
the  effect  of  additive  superposed  structures  can  be  seen,  an  effect 
caused  by  Interband  transitions.   As  expected  from  the  reflectivity, 

the  real  part  of  the  dielectric  constant  remains  negative  until  15  eV, 

■ 
when  the  reflectivity  begins  to  decrease.   The  imaginary  part  is  seen 

to  be  very  small  when  e  passes  through  zero  which  is  the  condition  for 

strong  plasma  resonance.   Figure  2.5  also  shows  the  conductivity 

CT  =  u)  €  /4tt  and  the  so-called  energy  loss  function  Im  e  .   The  latter 

function  represents  the  energy  lost  by  fast  electrons  traversing  a 

thin  film  of  a  metal  and  can  be  shown  to  be  very  sharply  peaked  about 

the  plasma  frequency  when  the  conditions  foi*  plasma  oscillations  are 

i  -1 

fulfilled.   The  peak  value  in  the  curve  for  Im  e   in  Figure  2.5  is 

calculated  at  15.2  eV  and  is  in  good  agreement  with  characteristic 

energy  loss  experiments. 

The  fact  that  in  aluminum  Interband  transitions  do  not  set  In 

until  about  1  eV  makes  possible  the  separation  of  the  free-electron 

contribution  from  the  Interbond  transitions.   It  is  therefore  possible 

to  fit  the  experimental  results  at  low  energies  to  Equation  (31). 

The  results  of  this  fitting  to  the  experimental  values  of  e  and  €_ 

allow  the  calculation  of  the  plasma  frequency  and  relaxation  time 

associated  with  the  conduction  band.   At  0.4  eV  the  calculated  values 

are:  htu„   =  12.7  eV  and  T  =  5.12  x  10~   sec. 
fa  c 
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Figure  2.5   The  spectral  dependence  of  the  real  and  Imaginary  parts 
of  the  dielectric  constant,  the  conductivity  a,   and  the 
energy  loss  function  Im  €~^  for  aliuninum.   (From  Ref- 
erence 14.  ) 
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It  Is  also  possible  to  compute  the  dc  conductivity.   For  the 
Drude  model  the  dc  conductivity  is  shown  to  be 


IJ0„   T 
Pa  c 

4n 


(34) 


17    -1 
Substitution  into  this  equation  yields  a    =  1. 55  x  10   sec   in  com- 

17    -1 
parison  to  the  experimental  value  3.18  x  10   sec   for  bulk  aluminum. 

■'■ 
This  rather  large  difference,  a  factor  of  about  2,   is  attributed  to 

the  fact  that  reflectivity  measurements  pertain  to  thin  films,  or 
alternatively  that  T  decreases  rapidly  from  its  dc  value  as  the  fre- 
quency is  Increased. 

The  value  of  the  optical  mass  is  also  determined  from  u) 

pa 

assuming  there  are  three  free  electrons  per  atom.   The  value  is 
m  =1.5.   Since  the  calculated  energy  bands  are  so  free-electron-like 
it  Is  clear  that  they  are  expected  to  lead  to  a  smaller  m  .   The  esti- 
mated optical  mass,  m   ,  predicted  by  band  theory  yield  a  value  close 

■  ■  ■ 
to  1. 15.   The  discrepancy  between  the  two  values  is  thought  to  arise 

from  the  interactions  neglected  in  the  individual  electron  model. 

(b) 
Since  e    =  0  in  the  region  near  zero  frequency  where  e  is 

largest,  the  separation  into  intra-  and  interband  parts  is  most  easily 

achieved  by 


calculating   e         from  e     (exp)   -  e         and  then  using 


Kramers-Kronlg  analysis   to   calculate   C 


The 


(b) 


results  for  the  interband  contribution  to  the  dielectric 
constant  are  shown  in  Figure  2.6.   These  curves  show  quite  clearly 
that  the  important  interband  transitions,  giving  rise  to  structure  in 
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Figure  2.6   The  real  and  imaginary  parts  of  e^  ^  (w)   for  aluminum. 
(From  Reference  14.) 
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the  reflectivity,  are  confined  to  a  small  energy  region  surrounding 
1.5  eV.   The  rapid  decrease  in  magnitude  of  both  dielectric  constants 
at  higher  energies  indicates  that  any  further  interband  transitions 
can  make  only  a  small  contribution. 


In 
convenient 


order  to  understand  the  Interband  effects  in  detail ,  it  is 
to  utilize  the  results  of  the  energy  band  calculations  for 


this  metal.  The  energy  bands  for  tho  various  symmetry  axes  are  shown 
in  Figure  2.7.  It  can  be  shown  [14]  from  Equation  (33)  that  the  con- 
tribution to  e    for  a  given  to  is  proportional  to  the  joint  density 

1(2  -• 

of  states,  times  an  average  of  \Pj>a'W      over  the  surface  tu-./ (k)  =co. 

Structure  Ir.  e  is  then  expected  to  occur  at  critical  points,  i.e.  , 
where  ^.^te'   =  0.   The  critical  points,  in  turn,  generally  occur  nt 
symmetry  points  where  the  gradients  of  E «(k)  and  E ./ (k)  may  vanish 
separately  due  to  symmetry. 

The  band  calculations  of  Figure  2.6  show  that  the  only  inter- 
band transitions  below  8  eV  at  symmetry  points  are  those  between  the 

two  occupied  states  W  'and  W  aiid  the  unoccupied  state  W  .   The  calcu- 
^      o  1 

lated  transition  energies  are  1.4  and  2.0  eV,  respectively.   It  can 
also  be  noted  that  interband  transitions  in  the  vicinity  of  D  also 
set  in  at  about  1.4  eV.   A  detailed  study  of  the  £(k)  in  this  region 
shows  that  a  critical  point  exists  on  the  [110]  axis  by  virtue  of  the 
fact  that  sl  and  Z  bands  below  and  above  the  Fermi  level,  respectively, 
are  parallel  there.   It  is  clear  then  that  the  regions  around  W  and  Z 
must  contribute  significantly  to  the  large  peak  in  e  . 
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Figure  2.7   The  calculated  energy  bands  of  aluminum  along  the 
symmetry  axes.   (From  Reference  14. ) 
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Band  calculations  predict  that  allowed  higher  energy  transi- 
tions at  symmetry  points  set  in  at  about  8  eV  and  above,  whereas  the 
experimental  e  exhibits  no  structure  in  this  range.   It  is  argued 

that,  for  the  case  of  aluminum,  the  peak  value  of  the  contribution 

-3/2 
from  a  critical  point  is  proportional  to  Eg     where  Eg  is  the  energy 

gap,   ITiis  fret  would  account  for  the  fact  that  no  strong  interband 

contributions  at  higher  energies  are  observed  for  this  metal. 

The  outstanding  feature  at  higher  energies  is  the  strong  and 
sharp  peak  in  the  energy  loss  function  Im  e~  at  15.2  eV.   This  energy 
is  generally  associated  with  the  plasma  resonance.   It  must  be  shown 
thot  this  value  is  consistent  with  the  result  huj    =  12.7  eV  obtained 
in  the  low^frequency  range  when  free-electron  effects  predominate. 

In  the  far  infrared  the  dielectric  constant  io  determined  by  the 

free-electron  expression  (31) ,  the  contribution  from  interband  processes 

being  relatively  unimportant.   Near  the  plosma  resonance,  u)  ,  at  about 

15  eV.  e„  is  very  small  so  that  U)  is  determined  quite  accurately  by 
2  p 

the  condition  that  e^  (lo  )  =  0.   This,  together  with  the  fact  that 
1  P 

lu  T  »  1  at  this  frequency,  implies  that 

2 


U) 


1  - 


pa 


uu 


^      1      P 


(35) 


and  hence  that 


U) 


pa 


""P   [l  +  ej''\u)  )  It 


']' 


40 


The  contribution  of  e^  ,  whose  value  is  about  -0.3,  is  therefore  quite 
important  in  determining  the  magnitude  of  U)  .  It  is  possible  to  deter- 
mine a  self-consistent  solution  of  Equation  (35)  ,  using  the  known  value 

of  CO   and  ef   as  given  in  Figure  2.6.   It  is  found  that  h  «)  =  15.2  eV, 
pa      1      •*         "  p         ' 

in  agreement  with  the  location  of  the  pemk  in  the  onergy  loss  function. 

The  fact  that  u)  is  greater  than  u)   results  from  the  negative  value  of 
P  pa 

ej   which  is  a  remnant  of  the  oscillator-like  dispersion  associated 
with  the  interband  transitions  near  1.5  eV. 

It  is  interesting  to  note  that  tho  dielectric  constant  can  also 

2    2 
be  approximated  by  its  general  asymptotic  form  e((u)  =  1  -  4TTne  /nu)  , 

where  m  is  the  free-electron  mass  and  n  corresponds  to  3  valence  elec- 
trons per  atom.   Thus,  the  plasma  frequency  can  be  written  in  the  form 

previously  used  in  connection  with  Equation  (21),  namely, 

2    i 
u)  =  (4Tine  /m)^.   Direct  substitution  into  the  preceding  expression 
P 

yields  h  ou  =15.8  eV,  in  good  agreement  with  the  earlier  result  of 
15.2  eV.   The  negative  value  of  e    is  therefore  seen  to  raise  the 

value  of  the  plasma  frequency  from  its  lower  value  (U   ,  as  determined 

2   A 
by  the  optical  effective  mass,  nearly  to  the  larger  value  (4Trne  /m)^ 

required  by  the  asymptotic  behavior. 

The  Noble  Metals.   Figure  2.8  shows  the  reflectivity  spectrum 
of  silver  to  energies  up  to  about  22  eV.   It  can  be  observed  that  at 
low  frequencies  the  reflectivity  is  near  100  percent,  already  pointed 
out  for  aluminum,  which  is  characteristic  of  the  free-electron  range. 
However,  near  4  eV  there  is  an  extremely  sharp  dip,  with  the  reflec- 
tivity dropping  to  less  than  1  percent  in  a  very  small  energy  range. 
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In  fact,  this  dip  Is  so  narrow  and  the  absorptivity  so  low  at  the 

bottom  of  the  dip  that  It  Is  possible  to  construct  a  narrow  bandpass 

I  ° 

filter,  using  a  thin  film  of  silver  (about  500  A)  that  will  transmit 

ultraviolet  light  in  the  vicinity  of  315  nm. 

At  energies  beyond  the  dip  the  structure  is  characteristic 

of  interbahd  processes.   Here  the  transitions  are  not  confined  to  a 

small  energy  region,  as  in  the  case  of  aluminum,  but  are  rather  widely 

dispersed,  as  occurs  with  most  metals.   The  reason  for  this  is  that  in 

the  case  of  silver,  and  the  noble  metals  in  general,  the  d  bands  are 

located  close  to  the  Fermi  level  and  have  oscillator  strengths  that 

■ 

extend  over  a  large  frequency  range. 

It  has  already  been  pointed  out  that  a  sharp  decrease  in  the 
reflectivity  is  usually  associated  with  a  plasma  resonance.   Figure  2.9, 
which  shows  the  real  and  imaginary  parts  of  the  dielectric  constant 
and  the  energy  loss  function  for  silver,  also  exhibits  a   sharp  peak 
in  the  energy  loss  function  near  4  eV.   In  this  range  c  goes  through 
zero  and  c  '  is  very  smalLas  required  for  the  existence  of  plasma 
oscillations.   The  data  shown  in  Figure  2.9  were  obtained  by  Kramers- 
Kronig  analysis  of  the  reflectivity  data.   For  this  analysis  the  reflec- 
tivity curye  was  extrapolated  beyond  22  eV  so  as  to  make  the  calculated 
values  of  e  and  e  agree  with  independent  measurements  of  Schulz  [16] 
and  Roberts  [17]  in  the  region  near  1  eV. 

It  can  also  be  observed  in  Figure  2.9  that  e     falls  off  rapidly 
in  the  free-electron  region,  as  expected  from  Equation  (31) ,  and 
approaches  a  value  very  close  to  zero  before  the  rise  due  to  the  onset 
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of  interband  transitions  follows.   This  fact  permits  a  relatively 
unambiguous  separation  of  e  into  its  "free"  and  "bound"  contributions. 
That  is,  to  a  good  approximation,  it  is  assumed  that  below  the  onset 

of  interband  processes  e    =0,  and  that  for  energies  above  the  crit- 

(f) 
ical  energy  e    =  0.   This  fact  may  be  used  in  conjunction  with 

Krnmers-Kronig  analysis  to  calculate  e.   .   Once  this  contribution  is 

separated,  e    is  obtained  from  e-   =  e  (exp)  -  e   .   The  results 

f 
of  this  analysis  is  shown  in  Figure  2.10.   The  curves  for  e  (cu)  are 

found  to  fit  accurately  the  description  of  the  Drude  region  given  by 

Equation  (31).   The  results  also  agree  quite  well  with  the  measurements 

of  Schulz  [16]  and  Roberts  [17],  giving  a  value  for  the  optical  mass 

m  =  1.03  and  a  relaxation  time  equal  to  that  obtained  from  dc  conduc- 
a 

tivity  measurements. 

In  connection  with  the  identification  of  plasma  oscillations, 

it  is  observed  in  Figure  2.9  that  in  the  absence  of  interband  transi- 

(f) 
tlons,  €  would  be  given  by  e    and  e  would  vanish  substantially  for 
1  1       ^ 

frequencies  greater  than  the  critical  frequency  to    The  plasma  reson- 
ance would  occur  where  e^  =  0,  i, e.  ,  co  =  u)   ,  the  value  calculated  for 

1    '  pa' 

the  density  of  free  electrons  having  mass  m  given  above.   Figure  2.10 

Indicates  that  h 'JU   =9.2  for  silver.   When  Interband  transitions 
pa 

contribute,  this  simple  picture  is  no  longer  valid,  since  the  contrlbu- 

.(b) 


tion 


of  e    is  such  as  to  cause  e  to  be  zero  in  a  region  where  e  is 

1  ^ 


also  small.   That  is,  the  increase  in  e  ,  due  to  the  onset  of  transitions 
from  the  d  band  to  the  Fermi  surface,  causes  e  to  vanish  at  an  energy 
Just  below  to  ,  the  threshold  for  these  transitions.   As  shown  in 
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Figure  2.8   Spectral  dependence  of  the  reflectivity  of  silver. 
(From  Reference  12.) 


A 


Figure  2.9   Spectral  dependence  of  the  real  and  imaginary  parts  of 

the  dielectric  constant  and  the  energy  loss  function 

2    2 
e^/(e  +  e  )  for  silver.   (From  Reference  12.) 
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Figure  2.10  Decomposition  of  the  experimental  value  of  e  for 


silver  into  free  and  bound  contributions 
(From  Reference  12.) 


(f)   ^   (b) 
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Figure  2.9,  this  occurs  when  e  happens  to  be  quite  small  and  is  thus 
responsible  for  the  sharp  and  well-localized  plasma  resonance  and  the 
sharp  peak  in  the  energy  loss  function  at  this  energy.   This  is  not 
a  free-electron  resonance  but  rather  a  hybrid  resonance  resulting  from 
the  cooperative  behavior  of  both  the  d  and  conduction  electrons. 

The  investigation  of  the  optical  properties  of  copper  follows 
very  much  along  the  lines  described  for  the  case  of  silver.   In  partic- 
ular, the  separation  of  the  dielectric  constant  into  "free"  and  "bound" 
contributions  is  carried  out,  using  the  same  technique  which  employs 
Kramers-Krpnig  analysis  and  the  experimentally  determined  values  of 
e^  and  c^  of  Schulz  [16]  and  Roberts  [17]  as  a  check  to  obtain  a  con- 
sistent analysis  of  the  data.   On  the  other  hand,  the  electronic  band 
structure  of  copper  has  been  worked  out  by  several  investigators  [18-20] 
along  the  symmetry  directions  in  k  space.   This  allows  the  interpreta- 
tion of  the  structure  in  the  optical  parameters  to  be  analyzed  in  more 
detail ,  particularly  In  regard  to  the  assignment  of  structure  to  partic- 
ular transitions. 

The  reflectivity  spectrum  of  copper  is  presented  in  Figure  2. 11 
for  er-^rgies  up  to  25  eV.   The  calculated  dielectric  constant  as  well 
as  the  energy  loss  function  obtained  from  the  reflectivity  is  shown 
in  Fi,?ure  2il2.   Among  the  outstanding  characteristics  of  the  data  for 
copper,  it  will  be  noted  that  interband  transitions  in  this  metal  set 
in  at  about  2  eV,  as  can  be  seen  from  the  sudden  rise  in  e  shown  in 
Figure  2.12.   This  in  turn  causes  the  reflectivity  of  copper  to  drop 
from  about  100  percent  typical  of  the  free-electron  ranges  producing 
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Figure  2.11   Spectral  dependence  of  the  reflectivity  of  copper. 
(From  Reference  12. ) 
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Figure  2.12   Spectral  dependence  of  the  real  and  imaginary  parts  of 

the  dielectric  constant  and  the  energy  loss  function 

2    2 
e„/(€  +  e  )  for  copper.   (From  Reference  12.) 
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a  sharp  shoulder  at  this  energy.   In  contrast  to  silver,  the  d  band 
in  copper  and  gold  lies  closer  to  the  Fermi  level  by  about  2  eV. 
In  the  latjter  metals  the  d  bands  are  located  about  2  to  2.5  eV  below 
the  Fermi  level.   Therefore,  the  onset  of  interband  processes  occur 
approximately  in  the  middle  of  the  visible  spectrum,  producing  a  drop 
in  the  reflectivity  which  is  responsible  for  the  yellowish  color  of 
copper  and  gold. 

Another  feature  of  the  optical  data  of  copper,  associated  with 

■ 

the  position  of  the  d  band,  is  the  fact  that  the  onset  of  interband 
tronsitions  takes  place  when  both  c  and  e  are  relatively  large  in 
magnitude  and  the  strength  of  these  transitions  is  not  large  enough 
to  force  6:  through  zero  at  a  smaller  energy.   This  point  is  illus- 
trated  in  Figure  2.13,  where  the  free  and  bound  contributions  to  the 
real  part  of  the  dielectric  constant  are  shown  separated.   The  free- 
electron  part  corresponds  to  m  =  1.42  and  a  relaxation  time  equal  to 


that  obtained  from  dc  conductivity  measurements,  namely, 
T„  =  3.5  X 


-14 
10    sec.   It  can  be  seen  that  with  the  free-electron 


contribution  only,  e  passes  through  zero  at  h  u)   =9.3  eV.   When 
the  interband  part  is  taken  into  account,  this  point  moves  to  about 
7.5  eV.   However,  the  conditions  for  strong  plasma  resonance  are  never 
fulfilled  in  copper  and  the  energy  loss  function,  shown  in  Figure  2.12, 
does  not  show  a  sharp  and  well-localized  plasma  resonance  like  the  one 
observed  in  silver. 

While  the  above  interpretation  of  the  optical  properties  of 
the  noble  metals  appears  to  be  essentially  correct,  the  original 
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Figure  2. 13  Decomposition  of  the  experimental  value  of  e     for  copper 
into  free  and  bound  contributions  e    and    e    . 
(From  Jteference  12.) 


49 


assignment  of  the  optical  transitions  made  by  Ehrenreich  and  co-workers 
(baaed  on  Segall's  band  calculations)  is  at  variance  with  results  of 

later  calculations.   Beaglehole  [21]  has  made  careful  measurements  on 

■ 
the, optical  parameters  of  copper  over  n  range  up  to  28  eV.   His  results 

are  in  good  agreement  with  those  of  Ehrenreich,  but  Beaglehole  suggests 

that  a  definite  assignment  to  particular  transitions  should  wait  for 

better  understanding  of  the  electronic  structure  of  copper.   On  the 

other  hand,  Mueller  and  Phillips  [22]  made  theoretical  calculations  of 

the  optical  spectrum  of  copper.   Their  results  appear  to  agree  well 

with  Burdick's  band  calculations  [19]  for  this  metal  and  furthermore 

predict  optical  transitions  in  agreement  with  suggestions  made  by 

Beaglehole.   It  should  also  be  pointed  out  that  the  energy  gaps  at 

symmetry  points  calculated  by  Segall  differ  from  those  calculated  by 

other  authors  [19,20]  by  as  much  as  20  percent.   It  would  appear  that 

Burdick's  band  calculations  on  copper  and  the  optical  density  of  states 

derived  from  his  work  fit  experimental  optical  data  a  great  deal  better 

than  do  the  calculations  of  Segall.   More  recently,  Fong  et  al.  [23] 

have  made  optical  studies  of  copper,  using  wavelength  modulation. 

They  have  also  made  band  calculations  [20]  that  agree  closely  with 

those  of  Burdick  [19]  and  the  interpretation  of  their  experimental  data 

is  consistent  with  the  band  structure. 

The  assignment  of  transitions  at  low  and  moderate  energies  in 

copper  appears  to  be  well  established  at  the  present  time.   Figure  2.14 

shows  the  electronic  band  structure  for  this  metal  as  obtained  by  Burdick. 

With  reference  to  the  structure  shown  in  Figures  2.11  and  2.12,  the 
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Figure  2.14  Calculated  band  structure  of  copper  along  symmetry 

axes.   Arrows,  indicate  interband  transitions  produc- 
ing structure  in  the  optical  constants.   (From 
Reference  19.) 
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sudden  shoulder  at  2-1  eV  is  associated  with  transitions  from  the  top 
of  the  d  band  ti  the  Fermi  surface  near  X  and  L,  whereas  the  structure 
at  about  4  eV  is  Identified  with  the  transitions  X_  to  X. /  and  L„/  to 
L-.   All  these  transitions  are  indicated  by  arrows  in  Figure  2.14. 

The  identification  of  optical  transitions  for  higher  energies 
is  more  difficult,  since  the  number  of  possible  transitions  is  greatly 
proliferated.   In  contrast  to  the  case  of  aluminum  and  most  semiconduc- 
tors, the  oscillator  strengths  in  the  noble  metals  appear  to  be  widely 
distributee!  in  energy.   Some  appreciation  of  this  fact  can  be  obtained 
from  consideration  of  the  sum  rule  yielding  n„ _,  the  effective  number 
of  electrons  per  atom  contributing  to  the  optical  properties  over 
a  given  frequency  range.   This  sum  rule  [10]  is  a  general  expression 
obtained  from  the  Kramers-Kronig  relations  which  can  be  written  as 


o  2    2 

u)  e^(uu)  dUJ  =  ^^  "  ^  (36) 

£•  m 


where  n  is  the  total  density  of  electrons.   If  the  core  states  in  the 
noble  metals  can  be  neglected  over  the  entire  far  vacuum  untraviolet 
range,  then  a  numerical  evaluation  of  n,  using  the  experimental  values 
of  £-  in  Equation  (36),  should  saturate  at  a  value  corresponding  to  11 
electrons  per  atom  (Is  plus  10  d  electrons)  and  the  sum  rule  is  said  to 
be  exhausted.   An  empirical  modification  of  Equation  (36)  can  be  made 
if  the  integration  is  carried  to  a  frequency  <o  rather  than  to  infinity. 
The  calculated  value  of  n  gives  therefore  a  measure  of  the  electrons 
taking  part  in  the  optical  processes  to  the  frequency  u)  and  the 
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rapidity  with  which  n   .  approaches  its  saturation  value  as  a  function 
of  energy  gives  some  information  concerning  the  distribution  of  oscil- 
lator strengths.   Figure  2.15  shows  n    as  a  function  of  energy 
obtained  from  numerical  integration  of  the  experimental  e  by  Ehrenreich 
et  al.  for  the  noble  metals  [12,15].   The  singularity  of  e^   at  zero 
frequency  causes  n    to  rise  very  rapidly  to  a  value  near  unity.   The 
deviation  from  unity  is  due  to  departures  of  the  effective  mass  value 
from  the  free-electron  mass.   Figure  2.15  shows  that  n    essentially 
saturates  at  the  free-electron  value  and  is  nearly  constant  until  the 
threshold  for  interband  transitions  is  reached.   The  fact  that  n  ._ 
Increases  very  slowly  and  is  for  copper  and  silver  still  far  removed 
from  the  expected  saturation  value  of  11  at  25  eV  indicates  that  the 
oscillator  strengths  for  transitions  involving  d  and  s  electrons  are 
distributed  over  a  rather  wide  energy  range. 
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Figure  2.15 


20   eV  U 


Effective  number  of  electrons  per  atom  as  a  function 
of  energy  obtained  from  numerical  integration  of  e 
for  silver,  coppefr  and  gold.   (Prom  Reference  9.) 


CHAPTER  III 
EXPERIMENTAL  PROCEDURE 

3.1  The  Apparatus 

The  differential  reflectometer  developed  in  this  work  is 
described  in  this  chapter.   The  instrument  fulfills  the  requirements 
set  forth  in  the  design;  namely,  the  measurement  of  small  differences 
in  reflectivity  between  two  specimens  at  room  temperature  in  the 
spectral  range  from  200  nm  to  700  nm.   The  discussion  can  be  sepa- 
rated conveniently  in  two  parts:   the  first  part  contains  the  modu- 
lating unit  and  various  optical  components;  the  second  part  is  formed 
by  circuits  and  instruments  arranged  for  the  reduction  of  the  data 
and  the  recording  of  the  changes  in  reflectivity  as  a  function  of  the 
incident  light.   The  technique  of  compositional  modulation  and  a 
first  version  of  the  apparatus  for  differential  reflectivity  mea- 
surements has  been  reported  [24]  in  a  preliminary  paper. 

Figure  3.1  shows  a  schematic  diagram  of  the  optical  part  of 
the  apparatus.   Light  that  comes  from  a  monochromator  passes  through 
lens  1  and  falls  onto  the  mirror  M  which  performs  small  oscillations 
at  the  line  frequency,  i.e.,  60  Hz.   As  a  consequence  the  light  re- 
flected from  this  mirror  illuminates  alternately  two  specimens  mounted 
side  jv  side  and  which  have  a  different  composition.   After  reflection 
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from  the  specimens  the  light  goes  through  lens  2  and  then  into  a  photo- 
multiplier  tube  (PMT). 

If  R^  and  R^  represent  the  reflectivities  of  each  of  the  speci- 
mens and  if  i^^  is  the  Intensity  of  light  at  a  particular  wavelength  \, 

then  the  light  reflected  from  the  specimens  will  have  intensities  I  R 

A  1 

^^^  ^\^2'      ^^®  reflectivities  of  the  sample  being  in  general  different, 
the  output  voltage  from  the  PMT  exhibits  small  steps  corresponding  to 
the  light  being  reflected  from  each  of  the  specimens.   This  is  indi- 
cated by  Cj  and  e  in  Figure  3. 1, 

The  geometry  of  the  optics  is  such  that  lens  1  images  the  exit 
silt  of  the  monochromator  onto  the  samples  so  that  the  light  spot  that 
scans  over  the  specimens  is  a   fine  line  parallel  to  the  boundory  sepa- 
rating the  alloys.   Lens  2,  on  the  other  hand,  is  placed  so  as  to  give 
a  magnified  image  of  the  vibrating  mirror  on  a  ground  plate  of  fused 
quartz  a  short  distance  In  front  of  the  face  of  the  photomultlplier 
tube.   This  plate  is  necessary  because  the  sensitivity  of  the  photo- 
cathode  is  not  uniform  over  its  entire  surface.   The  light  passing 
through  the  ground  plate  is  diffused  in  all  directions,  thereby  min- 
imizing the  effect  of  the  slight  motion  of  the  light  spot  that  results 
from  the  motion  of  the  oscillating  mirror. 

It  was  shown  In  preliminary  experiments  that  the  response  of 
the  PMT  to  the  incident  radiation  is  linear  over  a  wide  range  of 
intensities.   In  addition,  the  PMT  has  a  very  small  dark  current. 
Therefore  it  is  possible  to  express  its  output  voltage  e  as 


e  =  f(V,\)  A  I^   R,  (37) 
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where  A  Is  the  effective  area  of  the  exit  slit  and  f(V,\)  is  some 
unknown  function  of  the  wavelength  \,   which  represents  the  spectral 
sensitivity  of  the  detector,  and  the  voltage  V  applied  to  the  detector. 
It  therefore  follows  that  the  step  (e^  -  e  )  is  proportional  to  the 
difference  in  reflectivities  (R  -  R_)  between  the  two  alloys  and  can 
be  written  as 


«1  -  «2=  ^<V'^>  A  \   ^\   -V  • 


(38) 


A  high-pressure  xenon  lamp  was  used  ns  a  light  source,  which 
covers  the  spectrum  from  the  infrared  (1000  nm)  to  the  near  ultra- 
violet (200  nm)  with  an  output  intensity  at  200  nm  of  about  1  percent 
of  the  maximum  output  of  the  lamp. 

The  light  coming  from  the  light  source  Is  focused  into  the 
entrance  slit  of  a  plane  grating  monocheomator   that  covers  the  range 
from  200  nm  to  700  nm.   The  specifications  given  by  the  manufacturer  of 
this  grating  are  a  dispersion  of  7.4  nm/mm,  a  blaze  at  210  nm,  and 
a  typical  scattered  light  of  0.2  percent.   With  the  slits  opened  at 
about  0.3  mm  the  resulting  dispersion  at  the  exit  slit  is  approx- 
imately 2  nm  which  is  quite  adequate  considering  that  the  structure 
to  be  resolved  by  the  Instrument  is  a  great  deal  larger  than  the  spread 
in  wavelengths.   The  scattered  light,  however,  can  be  a  problem  at  the 


ultraviolet 


end  of  the  useful  spectrum,  because  here  the  output  from 


the  light  source  is  very  small.   This  may  be  appreciated  by  writing 


Model  613.   McPherson  Instrument  Corporation. 
Model  33-86-07.   Bausch  &  Lomb,  Inc. 
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for  the  light  Intensity  from  the  exit  slit 


i:  =  I  +  I,  ,  (39) 

A.     o     A. 


that  is,  the  output  is  the  sum  of  n  "pure"  monochromatic  component  I. 
plus  a  "white"  background  I  .   For  most  of  the  useful  range  the  light 
output  from  the  lamp  is  large  and  I   is  negligible.   Near  the  short- 
wavelength  end,  however,  I   is  quite  small  and  the  stray  light  becomes 
increasingly  disturbing.   Since  the  output  voltage  from  the  PMT  results 
from  both  terms  in  Equation  (39) ,  an  estimate  of  the  relative  magnitudes 
of  these  two  terms  can  be  obtained  by  removing  one  of  these  components 
from  the  light  beam  by  an  appropriate  filter.   In  the  ultraviolet, 
where  the  background  light  is  important,  this  can  conveniently  be  done 
by  inserting  n  glass  plate  into  the  light  beam.   This  plate  blocks  off 
the  ultraviolet,  allowing  most  of  the  background  to  pass  on  to  the  PWr. 
Therefore,  by  noting  the  readings  from  the  PMT  with  and  without  the 
glass  plate  for  various  wavelengths,  it  is  possible  to  measure  the 
fraction  of  the  PMT  output  which  is  due  to  stray  light.   Figure  3.2 
shows  the  results  of  measurements  carried  out  in  this  way.   It  can  be 
seen  that  at  200  nm  the  background  light  is  responsible  for  about 
15  percent  of  the  PMT  output  and  it  follows,  therefore,  that  results  at 
this  end  of  the  spectral  range  must  be  interpreted  with  caution. 

The  photomultiplier  is  of  the  S-20  "extended"  type  with  fused 
silica  entrance  window.   The  spectral  sensitivity  of  the  tube  covers 

EMI  9659. 
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the  range  between  about  900  nm  to  about  160  nm  which  matches  well  with 
the  monochroraator  range. 

Figure  3.3  shows  the  spectral  distribution  of  the  xenon  lamp 
multiplied  by  the  sensitivity  of  the  photomultiplier  tube  as  measured 
by  the  PMT  voltage  output.   It  can  be  seen  that  the  source  contains 
a  few  peaks  around  450-500  nm  which  is  characteristic  of  the  xenon 
emission  lines,  but  otherwise  it  is  a  fairly  continuous  light  source. 

The  vibrating  raix'ror  was  constructed  from  a  dc  moving  coll 
galvanometer  by  dismounting  the  face  and  bending  the  pointer  upward 
into  a  vertical  position.   A  small  front  face  aluminum  mirror  was  then 
attached  to  the  pointer.   When  a  small  ac  voltage  is  applied  to  the 
galvanometer,  it  vibrates  at  the  line  frequency.   The  amplitude  of  the 
vibrations  can  be  adjusted  by  changing  the  applied  voltage.   This 
amplitude  is  normally  set  so  as  to  produce  a  light  spot  on  the  speci- 
mens that  scans  over  a  span  of  approximately  3  mm.   In  addition,  the 
center  of  oscillations  can  be  varied  by  a  small  dc  bias  to  the  moving 
coil  which  allows  the  scanning  beam  to  be  centered  on  the  optical  axis. 
The  specimens  are  placed  on  a  moving  vertical  plate  that  can  be 
displaced  on  a  plane  parallel  to  the  reflecting  surfaces  by  means  of 
two  micrometer  screws.   This  is  done  in  order  to  select  the  areas  of 
the  specimens  that  are  to  be  examined  and  also  to  center  the  boundary 
between  the  specimens  on  the  scanning  beam. 

3.2  Data-Reducing  Circuit 

The  purpose  of  data-reducing  circuit  is  to  process  the  signal 
from  the  PMT  in  order  to  obtain  a  graph  of  the  changes  in  reflectivity 
as  a  function  of  wavelength. 
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Figure  3.3   Spectral  distribution  of  the  xenon  lamp  multiplied  by 
the  sensitivity  of  the  photomultlplier  tube. 
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The  output  e  from  the  PMT  can  be  split  Into  two  channels: 

a  dc  channel  with  a  signal  proportional  to  R  (where  R  is  the  average 

value  of  R^  and  R  )  and  an  ac  channel  proportional  to  AR.   The  ac 
12 

channel,  which  contains  the  modulation  in  the  reflectivity,  is  fed 
into  a  phase  sensitive  amplifier  (also  referred  to  as  lock-in  amplifier). 
This  amplifier  delivers  an  output  which  is  proportional  to  the  RMS  value 
of  the  fundamental  frequency  component  of  the  input  signal.   That  is, 
its  output  is  a  dc  signal  proportional  to  Ae,  or  equivalently  to  AR. 

This  signal  could  be  fed  directly  to  the  y-axis  of  a  recorder. 
However,  as  can  be  seen  from  Equation  (38),  this  graph  would  contain, 
in  addition  to  the  signal  AR  of  interest,  the  structure  in  the  light 
source  through  the  factor  I  and  also  the  particular  spectral  sensitiv- 
ity  of  the  detector.   Moreover,  fluctuations  in  the  light  source, 
which  are  always  present,  show  up  in  the  graph  as  background  noise. 

It  is  very  desirable,  therefore,  to  normalize  the  signal  by 
dividing  the  output  from  the  lock-in  amplifier  by  the  average  value 
e  =  (e^  +  e_)/2.   This  way  the  factor  f(V,\)  A  I   cancels  out. 
It  follows  that 

^  -  ^2  ^  "l  -"2 
e         R 

Problems  associated  with  the  instability  of  the  source  are  avoided  if 
the  normalization  is  performed  with  an  "on-line"  computer  element. 
However,  it  is  not  easy  to  construct  a  simple  electronic  analog  circuit 
to  carry  out  this  division.   The  large  dynamic  range  required  and  over- 
loading of  the  PMT  are  two  of  the  problems  that  arise. 
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One  possible  analog  system  [25]  is  obtained  by  feeding  the  sig- 
nals from  the  lock-in  amplifier  and  the  PMT  into  two  logarithmic  con- 
verters and  the  outputs  of  the ; converters  to  a  differential  amplifier 
so  as  to  obtain  log  AR/R. 

A  better  solution  to  the  normalization  hns  been  developed, 
using  operational  amplifiers  and  a  divider  component.   A  block  diagram 
of  the  circuit  is  shown  in  Figure  3.4.   The  PMT  output  is  fed  into  the 
lock-in  which  picks  up  and  amplifies  the  modulated  part  and  becomes 
the  numerator  channel.   The  PMT  output  is  also  fed  into  a  dc  amplifier 
which  includes  a  capacitor  in  its  feedback  loop.   The  signal  is  there- 
fore filtered  to  a  dc  signal  proportional  to  R  =  (R-  +  Rg)/^  becoming 
the  denominator  channel.   This  denominator  amplifier  is  wired  so  as  to 
satisfy  several  other  requirements:   it  presents  a  very  high  impedance 
to  the  PMT,  drawing  very  little  current  from  it;  it  changes  the  sign 
of  the  PMT  signal,  converting  it  into  a  positive  value  as  required  by 
the  divider  circuit;  it  also  permits  the  adjustment  of  the  voltage  gain 
in  order  to  operate  the  divider  in  its  optimum  range. 

The  attenuator  in  front  of  the  lock-in  amplifier  (see  Fig- 
ure 3.4)  is  necessary  because  of  the  relatively  high  level  signal  from 
the  PMT.   The  maximum  input  level  that  can  be  fed  into  the  lock-in 
amplifier  is  50  mV,  whereas  the  PMT  output  is  typically  a  few  volts. 
The  output  from  the  lock-in  could  be  fed  directly  into  the  numerator 
of  the  divider.   It  is,  however,  convenient  to  reduce  the  level  of  the 
lock-in  output  to  match  the  input  requirements  of  the  divider.   The  dc 
amplifier  shown  in  the  numerator  channel  acts  therefore  as  an  attenua- 
tor, the  gain  of  which  can  be  adjusted  by  a  potentiometer  in  the  feed- 
back loop. 


64 


a, 


3, 

•H 


65 


Figure  3.5  shows  the  circuit  diagrams  for  the  dc  amplifiers 
and  the  divider  circuit.   The  function  performed  by  each  amplifier 
is  also  indicated  next  to  each  circuit. 

The  divider  output  Is  proportional  to  the  ratio  of  the  two 
input  signals, 

R-  -  R., 


R 

and  this  is  connected  to  the  y-axis  of  an  x-y  recorder.   On  the  other 
hand,  a  potentiometer  is  attached  to  the  monochromator  In  a  way  that 
when  the  monochromator  scans  over  its  entire  wavelength  range,  a  small 
dc  voltage  proportional  to  \  Is  produced  to  drive  the  x-axis  of  the 
recorder.   The  recording  operation  Is  therefore  performed  in  an  auto- 
matic way. 

In  order  to  monitor  the  PMT  signal,  an  oscilloscope  Is  con- 
nected directly  across  the  PMT  output.   This  permits  readings  of  the 
voltage  outputs  and  a  direct  determination  of  the  percentage  varia- 
tion of  th^  reflectivity  of  the  samples. 

Since  the  output  from  the  divider  is  the  ratio  of  two 
processed  signals,  attention  must  be  paid  to  the  time  constants  intro- 
duced into  'each  signal.   The  filtering  capacitor  in  the  denominator 
channel  introdjces  a  time  constant  to  this  signal.   On  the  other  hand, 
the  lock-In  amplifier  in  the  numerator  channel  has  an  adjustable  time 
constant  which  can  be  varied  in  steps  from  0,1  m  sec  to  30  sec.   For 
a  proper  performance  of  the  Instrument  it  is  convenient  to  have  these 
two-time  constants  equal.   An  indication  of  any  mismatch  in  time 
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constants  can  be  obtained  by  suddenly  increasing  the  applivjd  voltage 
to  the  PMT  and  noting  the  motion  of  the  recording  pen.   If  the  time 
delay  introduced  into  each  channel  is  equal  for  both,  then  no  motion 
of  the  recording  pen  is  observed.   On  the  other  hand,  if,  say,  the  numer- 
ator channel  has  a  larger  time  constant,  the  pen  will  move  to  small 
values  upon  increasing  the  applied  voltage.   Therefore,  the  mismatch 
can  be  eliminated  by  adjusting  the  lock-in  time  constant  until  no 
motion  in  the  recording  pen  occurs  when  the  applied  voltage  is  suddenly 
changed. 

Perhaps  the  most  outstanding  advantage  of  the  divider  is  that 
the  applied  voltage  to  the  PMT  can  be  adjusted  during  recording.   This 
is  very  convenient  because  the  light  intensity  from  the  source  is  very 
small  in  the  very  important  ultraviolet  region  of  the  spectrum. 
In  fact,  the  output  from  the  inonochromator  tails  off  rapidly  near 
200  nm  (Figure  3.3).   By  increasing  the  applied  voltage  In  the  ultra- 
violet region  It  is  possible  to  keep  a  nearly  constant  output  from  the 
PMT,  thereby  lifting  the  signal  from  the  background  noise.   Moreover, 
the  divider  is  kept  within  a  small  dynamical  range  where  its  perfor- 
mance is  optimum.   This  permits  the  extension  of  tha  useful  range  of 
the  Instrument  to  about  200  nm  which  corresponds  to  almost  the  limit 
that  can  be  achieved  without  vacuum. 

If  the  applied  voltage  to  the  PMT  were  to  be  varied  so  as  to 
keep  the  output  precisely  constant,  using  a  servomechanism,  then  the 
denominator  would  be  constant  at  all  times  during  recording.   This 
means  that  the  divider  circuit  would  not  be  necessary.   This  is  in  fact 
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another  way  to  obtain  AR/R  [26].   This  is  not,  however,  a  very  simple 
method  and  its  design  requires  close  attention  to  transient  effects 
within  the  power  supply. 

3.3  Sample  Preparation 

The  alloys  are  prepared  from  high  purity  metals,  encapsulated 

—5 
in  Vycor  glass  in  a  vacuum  of  approximately  10   Torr  and  melted  in  an 

Induction  furnace.   After  coolins  the  specimens  are  cold  worked  and 
homogenized  for  about  two  days  at  a  temperature  somewhat  below  the 
solidus  line.   In  order  to  make  reflectivity  measurements  with  the 
instrument,  two  specimens  having  different  compositions  are  cold 
mounted  side  by  side  with  no  gap  In  between.   The  specimens  are  then 
mechanically  polished  with  standard  metallographic  techniques  down  to 
Ip,  diamond  compound.   Reflectivity  measurements  are  carried  out  within 
a  few  minutes  after  final  polishing. 

Since  mechanical  polishing  leaves  a  layer  of  cold  worked 
material  on  the  surface  of  the  specimens,  preliminary  measurements 
were  conducted  in  order  to  get  some  appreciation  of  the  effect  of 
mechanical  polishing  on  the  reflectivity.   Two  specimens  of  high  pur- 
ity copper  were  mounted  side  by  side;  one  was  heavily  cold  worked, 
while  the  other  was  annealed.   The  composite  was  mechanically  polished 
in  the  way  described  above  and  subsequently  electropolished  to  remove 
a  few  microns  of  material,  thereby  leaving  the  two  specimens  with  the 
same  surface  treatment  but  one  cold  worked  and  the  other  annealed. 
Examination  of  these  specimens  with  the  differential  reflectometer 
revealed  no  significant  difference  in  reflectivity  over  the  range  of 
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the  instrument.   Therefore,  mechanical  polishing  appears  to  be  a  satis- 
factory means  of  surface  preparation  at  least  in  the  range  from  200  nm 
to  700  nm.  There  are  indications,  however,  that  in  the  vacuum  ultra- 
violet surface  preparation  becomes  very  critical  [21], 

It  should  be  mentioned  at  this  point  that  the  differential 
reflectivity  measurements  are  made  on  two  samples  which  have  similar 
composition.   Since  these  samples  receive  the  same  polishing,  the 
error,  if  any,  due  to  surface  distortion  (and  also  due  to  oxidation) 
is  approximately  the  same.   Even  if  the  absolute  values  of  the  reflec- 
tivities are  changed  by  surface  effects,  the  difference  in  their 
reflectivities  Is  probably  not  affected  by  significant  amounts. 


CHAPTER  IV 


RESULTS  AND  DISCUSSION 


4.1  General  Remarks 

The  results  obtained  by  compositional  modulation  of  copper-based 
alloys  containing  additions  of  nickel,  aluminum  and  zinc  are  presented 
and  discussed  here.   The  alloys  of  copper  with  aluminum  and  zinc  have 
been  kept  within  the  cv-phase  region.   At  first  a  general  description 
on  the  main  features  of  a   "differential  spectrogram"  of  two  alloys 
containing  different  amounts  of  a  second  phase  X  will  be  given.   In 
Figure  4,1a  the  reflectivities  of  two  olloys  are  schematically  plotted 
versus  the  photon  energy  of  the  Impinging  light.   The  spectral  reflec- 
tivities of  the  alloys  show  a  marked  decrease  in  the  low  energy  region, 
as  copper  does.   In  this  schematic  diagram  it  is  assumed  that  the  onset 
of  this  decrease  (marked  A  and  C,  respectively)  is  shifted  in  alloy  2 
towards  higher  energies. 

The  differential  spectrogram  (M/R  versus  photon  energy, 
Figure  4. lb)  shows  maxima  and  minima  at  energies  where  the  difference 
between  the  two  reflectivities  is  largest,  i.e.,  at  places  where 
"structure"  in  the  spectral  reflectivity  is  found.   For  example, 
a  maximum  in  m/R   can  be  seen  at  that  energy,  where  the  slope  of  the 
spectral  reflectivity  is  largest.    In  the  following,  actual  results 
of  differential  spectrograms  are  given  and  discussed. 


Here  R  stands  for  the  average  of  R^^  and  H^ 

70 


71 


Photon  energy 


Figure  4.1   (a)  Schematic  representation  of  the  spectral  reflectivities 
of  two  alloys  (marked  Alloy  1  and  Alloy  2).   (b)  Differ- 
ential spectrogram  for  the  two  alloys  referred  to  in  (a). 


72 


4.2  Copper-Zinc  Alloys 

Reflectivity  measurements  were  carried  out  in  a-copper-zinc 
alloys  with  the  aim  to  Investigate  the  changes  in  its  band  structure 
as  a  function  of  composition,  in  particular  the  behavior  of  the  main 
absorption  edge  around  2.1  eV  and  absorptions  near  4.0  and  4.4  eV. 
This  alloy  system  is  of  particular  interest  because  parts  of  the  band 
structure  for  o-brass  are  known  [27]  (Figure  4.2)  which  provides  a  baisis 
for  comparison  with  our  experimental  data.   For  calculation,  the  Kohn- 
Rostoker  method  in  combination  with  a  "pseudoper iodic"  potential  was 
used.   The  bands  of  a  copper-zinc  alloy  containing  30  at.  %  Zn  have 
been  calculated  [27]  at  high  symmetry  points  of  the  Brillouin  zone  with 
the  assumption  that  the  coppered  bands  ore  unchanged  by  alloying.   If 
one  compares  the  bands  of  pure  copper  as  calculated  by  Segall  [18]  with 
the  alloys  bands  (Figure  4.2),  one  can  observe  a  similar  shape.   The 
Fermi  level,  however,  has  changed  relative  to  the  d  bands.   This  can 
be  understood  by  considering  the  fact  that  zinc  has  essentially  one 
extra  electron  compared  to  copper.   This  has  the  effect  of  raising  the 
Fermi  level  by  filling  empty  stctes  (rigid  band  behavior).   A  predom- 
inant effect  due  to  alloying  is  therefore  an  increase  of  the  energy  gap 
between  the  top  of  the  d  bands  and  the  Fermi  surface.   Of  significance 
also  is  the  relatively  large  displacements  of  the  L  and  X  /  levels 
because  of  the  shift  of  the  s  and  d  bands  to  lower  energies.   This 
results  in  reduced  L_/  -L  and  X  -X  /  energy  gaps. 

Figure  4,3  shows  four  differential  spectrograms  (AR/R  vs  photon 
energy)  for  Cu-Zn  alloys  varying  from  a  medium  composition  of  2.25  at.  % 
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Figure  4,2       Comparison  of  energy  bands  of  Cu-30   at.  %    Zn   (dashed 
curves)   with  those  of   copper   (solid  curves).       (From 
Reference  27.  )  ] 
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Zn  to  11  at.  %  Zn.   The  energy  where  M/R  starts  to  increase  (marked 

T  in  Figure  4.3)  is  2.1  eV.   We  interpret  vhis  energy,  as  the  threshold 

energy  for  optical  Interband  transitions.- 

The  maximum  marked  A  in  the  differential  spectrograms  is  shifted 

to  larger  energies  with  increasing  concentrations  of  zinc.   We  ascribe 

this  peak  to  the  A_  -  A,  transition   (Figure  4,2).   According  to  the 
5    1 

considerations  mentioned  above,  this  energy  gap  increases  with  increas- 
ing zinc  concentration  due  to  a  raise  of  the  Fermi  level.   For  a-brass 
with  25  at.  7f  Zn  a  total  displacement  of  point  A  of  0.5  eV  has  been 
measured.   The  calculated  increase  of  this  energy  gap  was  found  to  be 
0.48  oV  [27]  which  is  in  good  agreement  with  our  optical  measurements. 

The  minimum  B  and  the  maximum  C  in  the  differential  spectro- 
grams of  Figure  4.3  are  shifted  to  lower  energies  with  increasing  zinc 
concentration  as  the  L^'  -  L-  and  the  X   -  X  /  energy  gaps  (see  above). 
We  associate  therefore  peaks  D  and  C  with  these  critical  interband 
transitions. 

It  should  be  emphasized  at  this  point  that  there  is  no  way  of 
knowing  exactly  that  a  certain  maximum  or  minimum  in  a  differential 
spectrogram  is  caused  by  a  very  specific  interband  transition.   It  is 
known,  however,  that  optical  interband  transitions  cause  structure 
in  the  function  c  = f (E)  and  therefore  most  probably  cause  structure 
in  R=f(E)  (see  Chapter  II).   Any  structure  in  the  spectral  reflectiv- 
ity curve  shows  up  in  a  differential  spectrogram  AR/R  =  f(E)  as  an 
extremum  point.   This  may  be  seen  by  considering  the  equation 


* 
And  all  the  transitions  from  the  5th  band  to  the  Fermi  surface. 


Cu-Zn 
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Figure  4.3   M/R  vs  wavelength  for  various  copper-zinc  alloys.   The 

first  number  on  the! lower  left  corner  of  each  graph  gives 
the  zinc: concentration  of  specimen  1  in  at.  %,   while  the 
second  number  gives  ,the  corresponding  information  for 
specimen  2.  Positive  values  of  ^/R  indicate  that  speci- 
men 2  has  higher  reflectivity. 
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de  (to)   dE   de  ((ju) 


where  e^   is  the  imaginary  part  of  the  complex  dielectric  constant 

(which  is  directly  related  to  R)  and  E   is  the  critical  energy  for 

cr 

optical  interband  transitions.   The  compositional  difference  between 
the  two  samples,  AC,  is  assumed  to  be  small  and  constant.   The  first 
factor  on  the  right-hand  side  of  Equation  (40)  is  independent  of  lu  and 
is  a  constant  for  a  particular  material  and  transition.   With  increas- 
ing photon  energy  hu)  ,  the:  second  term  is  expected  to  peak  at  E   , 
since  e^  usually  exhibits  a  singularity  at  this  energy  [28]. 

4.3  Copper-Aluminum  Alloys 

Another  interesting  system  to  study  changes  of  the  band  struc- 
ture by  alloying  is  represented  by  the  copper-aluminum  alloys.   The 
band  structure  of  aluminum  is  well  understood  and  Cu-Al  alloys  can  be 
considered  to  be  a  prototype  of  a  noble  metal  alloyed  with  a  metal  of 
nearly  free-electron-like  band  structure. 

Figure  4.4  shows  data  for  the  change  in  reflectivity  upon 
adding  1  at.  %  Al  to  Cu.   As  in  the  previous  case,  one  can  distinguish 
the  threshold  energy  T  and  various  peaks  designated  A  through  C. 
These  peaks  broaden  as  aluminum  is  alloyed  to  copper. 

In  Figure  4.5  the  shifts  of  the  aforementioned  peaks  are  pre- 
sented as  a  function  of  composition.   It  can  be  seen  that  similar  to 
the  case  of  or-Cu-Zn  alloys,  peak  A  (which  we  ascribed  to  the  A  -  A 
transition)  is  shifted  to  higher  energies  when  aluminum  is  added  to 
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Figure  4.4   Spectral  change  In  reflectivity  upon  adding  1  at.  %  Al 
to  Cu.   Positive  values  of  AR/R  indicate  that  the  alloy 
has  higher  reflectivity  than  pure  copper. 
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Figiire  4.5   Shifts  in  the  peaks  marked  A,  B  and  C  in  Figure  4.4  as 
a  function  of  concentration  of  aluminum. 


79 


copper,  whereas  peaks  B  and  C  are  shifted  to  lower  energies.   The 
interpretations  of  these  shifts  have  to  follow  a  similar  line  as 
described  in  the  previous  section. 

4.4   Copper-Nickel  Alloys 

Often  regarded  as  characteristic  of  the  many  alloy  systems 
between  a  noble  metal  and  a  transition  metal  are  the  Cu-Ni  alloys. 
They  are  ferromagnetic  over  more  than  one-half  of  the  composition 
range  and  are  substitutional  solid  solutions  over  the  entire  range. 
These  considerations,  and  the  fact  that  copper  and  nickel  are  among 
the  best  understood  of  the  noble  and  transition  metals,  has  stimulated 
research  of  the  Cu-Ni  system.   Early  measurements  on  Cu-Ni  alloys  led 
Mott  to  propose  the  rigid-band  model.   This  model  assumes  that  there 
is  one  electronic  density-of-states  function  which  is  the  same  for 
copper,  nickel,  and  Cu-Ni  alloys,  with  this  density  of  states  filled  to 
an  energy  level  determined  by  the  ratio  of  electrons  to  atoms.   Many 
subsequent  measurements,  however,  have  questioned  the  validity  of  the 
rigid-band  model  for  Cu-Ni  alloys  and  alternative  models  have  been 
suggested.   In  particular,  the  virtual-bound-state  model,  developed  by 
Frledel  [29]  and  Anderson  [30]  appears  to  be  applicable  to  the  Cu-rich 
Cu-Ni  alloys. 

These  two  models  predict  quite  different  behavior  for  the 
density  of  states  in  Cu-Ni  alloys.   This  is  illustrated  in  Figure  4.6. 
The  drawing  on  the  left  illustrates  schematically  the  density  of  states 
of  pure  copper  which  is  characterized  by  a  high  density  of  d  states  for 
energies  greater  than  2  eV  below  the  Fermi  level  (E_) ,  and  a  relatively 
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Figure  4.6   Illustration  of  the  density  of  states  in  copper  and 
copper-nickel  alloys,  showing  the  behavior  expected 
from  rigid-band  and  virtual -bound-state  models  at 
15  at.  %  Ni  in  Cu. 
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lower  density  of  s  and  p  states  between  the  d  bend  nnd  E  .   The  rigid- 

F 

band  model  for  Cu-Nl  alloys  assumes  that  when  nickel  is  added  to 
copper,  electrons  from  the  s-  and  p-derived  states  fill  the  nickel  d 
shell,  thereby  decreasing  the  separation  between  the  Fermi  level  and 
d  band.   For  15  ot.  %  Ni  in  copper,  the  alloy  density  of  states  would 
be  as  shown  schematically  in  the  central  drawing  of  Figure  4.6,  the 

separation  between  E-  and  the  upper  edge  of  the  d  bands  would  be 

F 

decreased  to  about  1,4  eV,  the  total  density  of  states  remaining  essen- 
tially unchanged. 

The  virtual-bound-state  model,  on  the  other  hand,  assumes  that 
the  impurity  3d  orbitals  of  nickel  form  levels,  highly  localized  around 
the  impuri».y  atoms,  that  are  broadened  in  energy  around  E. ,  the  center 
of  the  state,  through  a  resonant  scattering  interaction  with  the  nearly 
free-electron-like  band  of  the  host  metal  copper.   The  resultant  density 
of  states  expected  is  shown  schematically  in  the  right  of  Figure  4.6  for 
15  at.  %  Ni  in  copper.   Copper  and  nickel  d  electrons  would  form  essen- 
tially independent  levels,  the  energy  separation  between  copper  d  states 
and  E  would  be  unchanged,  and  the  nickel  d  electrons  would  cause  an 
increase  in  the  alloy  density  of  states  between  the  copper  d  band  and 
Fermi  surface. 

To  test  the  validity  of  one  or  the  other  of  these  two  models, 
the  following  experiments  were  undertaken.   Reflectivity  measurements 
with  the  differential  reflectometer  were  carried  out  in  copper-nickel 
alloys  in  the  composition  range  from  0.5  at.  %  Ni  to  about  15  at.  %  Ni. 
In  Figure  4,7  six  differential  spectrograms  for  various  copper-nickel 
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alloys  are  shown.   It  can  be  seen  that  the  edge  at  2.1  eV  (590  nm) 
in  pure  copper  does  not  move  to  lower  energies  when  nickel  Is  added 
to  copper  as  would  be  expected  on  the  basis  of  the  rigid-band  model. 

The  virtual-bound-stnte  model,  on  the  other  hand,  predicts  that 
this  edge  is  unaffected  by  alloying.   Moreover,  this  model  predicts 
a  drop  in  the  reflectivity  at  the  red  end  of  the  spectrum  from  the 
near  100  percent  value  in  pure  fsopper.   The  reflectivity  measurements 
(see  Figure  4,7)  show  that  this  dj  )p   does  take  place;  from  the  graph 
marked  0/10  it  can  be  seen  that  the  reflectivity  at  low  energies 
decreases  by  about  20  percent  with  the  addition  of  10  at.  %  nickel 
to  copper.   This  is  expected  because  the  d  states  from  nickel,  local- 
ized in  the  gap  between  the  copper  d  band  and  the  Fermi  surface,  provide 
a  higher  density  of  states  spread  in  the  energy  gap,  allowing  transi- 
tions to  take  place  between  these  bands  and  the  Fermi  surface,  thereby 
absorbing  energy  from  the  incident  radiation.   Theoretical  calculations 
made  by  Ehrenreich  and  co-workers  [31]  on  the  Cu-Ni  system  and  photo- 
emission  studies  of  this  alloy  made  by  Spicer  [32]  give  further  sup- 
port to  the  applicability  of  the  virtual-bound-state  model  to  copper- 
rich  Cu-Ni  alloys. 


The  wavelength  can  be  converted  to  energy  of  the  radiation 
with  the  expression  E  =  1240/X.   With  X  in  nm,  the  formula  gives 
E  in  electron  volts  (eV). 
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4.  5  Effect  of  Inhomogeneltles 

In  order  to  determine  the  amount  of  compositional  inhomogeneity 
of  the  samples  and  assess  Its  effect  on  the  optical  properties,  tests 
were  conducted  on  Cu-Nl  alloys  using  the  microprobe.   Two  specimens 
of  the  same  average  composition,  namely,  Cu-15  at.  %  Ni,  were  mounted 
side  by  side  in  the  manner  used  for  differential  reflectivity  measure- 
ments.  One  specimen  was  as  cast  and  the  other  cold  worked  and  homogen- 
ized at  about  1000°C  for  24  hours.   The  specimens  were  then  carefully 
polished.   Microprrbe  analysis  (wit   n  spatial  resolution  of  about  2^,) 
showed  that  there  was  a  mensurable  difference  in  composition  across  the 
surface  of  the  as  cast  sample.   The  homogenized  sample,  on  the  other 
hand,  showed  a  very  uniform  composition.   Tlie  chart  recording  from  the 
as  cast  sample  is  presented  in  Figure  4.8,  showing  a  compositional 
variation  between  79  at.  %  Cu  and  89  at.  %Cu.   Figure  4.9  shows  the 
corresponding  data  for  the  homogenized  sample. 

Of  particular  interest  was  the  determination  of  any  correla- 
tion between  the  microprobe  analysis  and  reflectivity  measurements. 
The  same  specimens  were  therefore  examined  with  the  differential 
reflectometer.   The  measurements  showed  a  difference  in  reflectivity 
between  1  and  2  percent.   The  results  are  shown  in  Figure  4.10  where  the 
higher  reflectivity  corresponds  to  the  as  cast  specimen. 

Having  these  results  in  mind,  particular  care  was  taken  in  the 
aforementioned  experiments  to  properly  homogenize  the  various  alloys 
before  measurements. 
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Figure  4.8   Compositional  variation  of  as  cast  Cu  -  15  at.  %  Ni 
alloy  as  revealed  by  microprobe  examination. 
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Figure  4,9   Compositional  variation  of  homogenized  Cu  -  15  at.  %  Ni 
alloy  as  revealed  by  microprobe  examination. 
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Figure  4.10   Spectral  change  in  reflectivity  between  as  cast  and 
homogenized  Cu  -  15  at.  %  Ni  alloys. 
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